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Introducing the Math Pac 


The Math Pac is a set of powerful tools for solving a wide range of mathematical, scientific, and en- 
gineering problems. These tools are provided in the convenient and flexible form of BASIC keywords. 
Once the math module is plugged into your HP-71 Computer, these keywords are instantly available: 
no program to load, no waiting. You can use these keywords in any program as often as needed; you 
avoid the restrictions that would apply to program calls and save the memory that subroutines would 
require. 


The Math Pac adds the following capabilities to your HP-71. 


© Complex variables and arrays. 


© Advanced real- and complex-valued functions. 


* Real and complex array operations. 


* Solutions to systems of equations. 


Roots of polynomial equations and user-defined functions. 


Numerical integration. 


© Finite Fourier transform. 
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How To Use This Manual 


This manual assumes that you are generally familiar with the operation of your HP-71 Computer, 
especially how to create, edit, store, and run programs. You should also understand the mathematical 
basis for the operations you will be performing. Because the keywords in the Math Pac cover such a 
wide range of mathematical subjects, we cannot provide much tutorial information on the mathematical 
concepts involved. 


The keywords in the Math Pac are independent of one another, so you may deal with only the 
keywords that specifically interest you. Each section in this manual contains information on keywords 
of a particular mathematical type—complex functions and operations, array arithmetic, and so on. All 
keywords described after section 5 (except FHF COT and INTEGRAL) use arrays in their operation. 
For an introduction to arrays, as used with the HP-71, read sections 3 and 14 of the HP-71 Owner's 
Manual. 


Variable Declarations 


The examples and programs in the Math Pac assume all variables are simple real unless otherwise 
declared. If an ERF;0(ats Tupe occurs as you execute an example or program, declare as REAL any 
variable not otherwise declared and continue operation. 


Array Types 


The Math Pac refers to two types of arrays, vectors and matrices. As used in this manual, the term 
vector identifies a singly-subscripted array, and matrix identifies a doubly-subscripted array. A sub- 
script must be a real numeric expression. At run time, a subscript expression is rounded to an integer. 
The value of this integer must be in the range [0,65535] (GF TION BASE &) or [1,65535] 
(OF TIOH BASE 1). Of course, in virtually all cases, available memory will determine the largest 
subscript you can use. 





An array can be one of five data types: F. ORT, INTEGER, COMPLEX, or COMPLE# 5 
(refer to section 3 for a description of COMPLEX and COMPLE SHORT). Math Pac MAT statements 
will not change the declared type of an array; for example, when the values from a REAL array are 
assigned to a SHORT or INTEGEF array, the values are rounded as they are stored into that array. 
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Array Redimensioning 


Some Math Pac keywords allow you to optionally redimension an array. This is called explicit 
redimensioning. Other keywords automatically redimension result arrays, if possible, to accomodate the 
number of elements generated by the keyword’s action. This is called implicit redimensioning. The kind 
of array redimensioning performed by a keyword, explicit or implicit, is stated in each keyword’s 
description. 


Explicit redimensioning occurs when an array’s size and subscript count is changed according to the 
number and value of new subscripts supplied by you. For example, if A is a 3 x 4 FEAL type matrix, 
then the HP-71 statement REAL A‘ 32> explicitly redimensions A to be a 3 dimensional vector. Note 
that explicit redimensioning allows arrays to be changed from vectors to matrices and vice-versa. Ex- 
plicit redimensioning also re-evaluates OP TIGH EASE; that is, resets the lower bound of an array’s 
subscripts if the OF TION EASE setting has changed. 





Implicit redimensioning occurs only in Math Pac operations of the form 
MAT result array = operation (operand array(s)). 


Implicit redimensioning only changes an array’s size. It does not allow changes between vectors and 
matrices, nor does it re-evaluate OF TIOH BASE. 


Keyword Description 


Within each section you will find a description of each keyword name, function, syntax, and operation 
in the following format. 


KEYWORD NAME Function That the Keyword Performs 





Syntax 


Legal data types and numeric values for use with this keyword. 








| Description of the values returned by this keyword and the details of the keyword’s operation. 


Keyword Name. This is the way the keyword will be referenced elsewhere in the manual. It is usually 
a mnemonic of the function that the keyword performs. In most cases the name must be embedded in a 
longer statement that includes arguments, parentheses, and so on; the name by itself usually isn’t an 
acceptable BASIC statement. 


Several keywords have names that are identical to names of keywords already present in your HP-71— 
like DISF, +, and #. The syntax in which such a name is embedded indicates which operation to 
perform. All operations available to you in the HP-71 itself are still available, unaffected by the pres- 
ence of the Math Pac. 
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Syntax. This is a description of the acceptable BASIC statements in which the keyword’s name can 
be embedded. The following conventions are used throughout the manual in describing the syntax of a 
keyword. 


Typographical Item Interpretation 


OOT MATRIX Words in dot matrix (like COMPLE) can be entered in lowercase or upper- 
case letters. The examples in this manual show statements, functions, and op- 
erators entered in UPPERCASE. 


italic Items in italics are the variables or parameters you supply, such as X in the 
SIMH¢X2> statement. 

bold Variables in bold type represent arrays. 

0) Square brackets enclose optional items. For instance, MAT A=ION[tX, Yo] 


indicates the redimensioning subscripts X and Y are optional. 
stacked items When items are placed one above the other, one and only one must be chosen. 


An ellipsis indicates that the optional items within the brackets can be re- 
peated. For instance, MAT INFUT A[, Bl]... indicates that MAT INPUT re- 
quires at least one array variable, and may accept several, with the array 
variables separated by commas. 


Legal Data Types and Numeric Values. This information, in the same box as the syntax, describes 
the types and ranges of arguments for the keyword that are acceptable to the Math Pac. Use this 
information to avoid generating errors and to isolate the cause of those that do occur. This is not a 
mathematical definition of the domain of the function that the keyword computes. 


Values Returned and Details of Operation. This information, in the box just below the syntax 
box, describes how the keyword works, tells what values the keyword returns, states whether array 
redimensioning (if any) is explicit or implicit, and states whether or not the keyword is usable in CALC 
mode. 


Examples 


Included in each section are a number of examples illustrating the use of the keywords in the section. 
To try an example yourself, type in the statements given in the Input/Result column using either 
upper- or lowercase, ending each line with with an (END LINE]. After you complete a line, the display of 
your HP-71 should look like the display shown in the Input/Result column following the line—pro- 
vided that you have set your HP-71 operating conditions as indicated below. 
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e All operating conditions should be set as listed in the reference manual in the Systems 
Characteristics Section under the topic Reset Conditions, except for those whose settings follow. 


e Set line width to 22 by entering WIDTH 22 [ENDLINE). 


e Set DELAY so that each display in a sequence of displays, often produced by a single statement, 
will remain visible long enough to be read and understood. The [EL fi’ statement is described in 
The HP-71 Reference Manual and section 1 of the HP-71 Owner’s Manual. In each you'll find 
descriptions of how you can control the length of time each display remains visible. For the display 
of array elements, you may find a DELAY & setting useful. This causes each display to remain 


until any key, such as [END LINE], is pressed. 


Additional Information 


Some sections in the Math Pac include additional information to help you make effective use of the 
more sophisticated operations. If you would like still more information, you can refer to the HP-15C 
Advanced Functions Handbook. Although the Math Pac differs from the HP-15C Advanced Pro- 
grammable Scientific Calculator in its operation and capabilities, much of the information in the 
HP-15C Advanced Functions Handbook applies to the Math Pac. Such information includes techniques 
to increase the effectiveness of equation-solving algorithms, integration algorithms, matrix operations, 
system solutions, and accuracy of numerical calculations. 


Section 1 
Installing and Removing the Module 


The Math Pac module can be plugged into any of the four ports on the front edge of the computer. 


CAUTIONS 


| Be sure to turn off the HP-71 (press [f] [OFF]) before installing or removing the module. 
| 
« If you have removed a module to make a port available for the math module, before installing the 





math module, turn the computer on and then off to reset internal pointers. 


* Do not place fingers, tools, or other objects into any of the ports. Such actions could result in minor 
electrical shock hazard and interference with pacemaker devices worn by some persons. Damage 
to port contacts and internal circuitry could also result. 


| 
* If a module jams when inserted into a port, it may be upside down. Attempting to force it further | 
may result in damage to the computer or the module. 
| 
| 
J 


| « Handle the plug-in modules very carefully while they are out of the computer. Do not insert any 
objects in the module connecter socket. Always keep a blank module in the computer's port when 
a module is not installed. Failure to observe these cautions may result in damage to the module or 
the computer. 


To insert the Math Pac module, orient it so that the label is 
right-side up, hold the computer with the keyboard facing up, and 
push in the module until it snaps into place. During this opera- 
tion be sure to observe the precautions described above. 





To remove the module, use your fingernails to grasp the lip on the bottom of the front edge of the 
module and pull the module straight out of the port. Install a blank module in the port to protect the 
contacts inside. 
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Section 2 
Base Conversions 


Binary, Octal, and Hexadecimal Representations 


The operations in this section allow your HP-71 to recognize and manipulate numbers expressed in 
number systems other than decimal (base 10). 


Because the HP-71 assumes that any real number stored in a numeric variable or entered from the 
keyboard is a decimal number, you must enter and store every non-decimal number as a character 
string. In particular, if you store the number in a variable, the variable’s name must end with “$”; if 
you enter the number from the keyboard, it must be enclosed in quotes. 


In the tables below, S$ will represent a binary, octal, or hexadecimal string or string expression. 


« A binary string consists entirely of 0’s and 1’s, and represents a number in the base 2 number 
system. A binary string expression is a string expression whose value is a binary string. 

e An octal string consists entirely of 0’s, 1’s, ..., 6's, and 7’s, and represents a number in the base 8 
number system. An octal string expression is a string expression whose value is an octal string. 


e A hexadecimal string consists of 0’s, ..., 9's, A’s, ..., and F’s (the letters may be either uppercase or 
lowercase), and this string represents a number in the base 16 number system. A hexadecimal string 
expression is a string expression whose value is a hexadecimal string. 


Base Conversion Functions 


BVAL Binary, Octal, or Hexadecimal to Decimal Conversion 


| pute S$,N3 
| 


where S$ is a binary string expression whose value is not greater than 
| 1110100011010100101001010000111111111111 (binary), and N is a numeric expression whose 
rounded integer value is 2; 


or S$ is an octal string expression whose value is not greater than 16432451207777 (octal), and N is a 
numeric expression whose rounded integer value is 8; 


or S$ is a hexadecimal string expression whose value is not greater than E8D4A50FFF (hexadecimal), 
and N is a numeric expression whose rounded integer value is 16. 
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BVAL (continued) 

r : = 

| Converts a string expression S$ representing a number expressed in base N into the equivalent decimal 
number. The value of the decimal equivalent can’t exceed 999,999,999,999 (decimal). 


Not usable in CALC mode. 


BSTR$ Decimal to Binary, Octal, or Hexadecimal Conversion 


BSTR#¢X,N23 


where X is a numeric expression, 0 < X < 999,999,999,999.5, and N is a numeric expression whose 
| rounded integer value is 2, 8, or 16. 


| Converts the rounded integer value of X (decimal) into the equivalent base N string. 
When N = 16, returns uppercase A, ..., F. 
Not usable in CALC mode. 











Examples 


dia 
VAL C "LBL", 


Se ig 


B4="1111" (ENDLINE] ee 


(BS, END LINE 





(END LINE] LINE 








i5 The decimal value of the binary string “1111.” 





+ LEND LINE 


The decimal value of the binary string 
“Tita 


END LINE 


i it 














nversions 17 





* (END LINE 


ile a arene oer 


BSTR#¢BVALC"AFICS",169,2> 
END LINE 




















The binary representation of AF1C8 
(hexadecimal). 


The octal sum of 14772 (octal) and 570 (octal). 





Additional Information 


Three considerations determined the range of acceptable parameters for the base conversion keywords. 
« The keywords give the exact answer for any integer in the range of acceptable parameters. 


* The keywords are inverses of one another, so that composition in either direction is the identity 
transformation for integers. 


The integers from 0 through 999,999,999,999 form the largest block of consecutive non-negative 
integers that the HP-71 can cisplay in integer format. 


Section 3 
Complex Variables 


Complex Data Types 


The operations in this section allow your HP-71 to declare, recognize and manipulate complex num- 
bers. These operations include: 
© Declaration of complex variables and arrays using COMPLE and COMPLE SHORT statements. 
° Extension of HP-71 variable assignment and the RES function to the complex case. 
° Extension of HP-71 IMAGE format strings to include complex fields 


* Conversion of real numbers to complex. 


Declaring Complex Variables 
COMPLEX Complex Variable Creation with 12-Digit Precision 


COMPLE variable list 


where the syntax is the same as that used for REAL, SHORT, and IHTEGER keywords. That is, each 
variable in the variable list has the form numeric variable [« dim 1[, dim 2] »], and dim 1 and dim 2 are 
real numeric expressions. 


Not usable in CALC mode. 





COMPLEX SHORT Complex Variable Creation with 5-Digit Precision 





COMPLEX SHORT variable list 





where the syntax is the same as that used for REAL, SHORT, and INTEGER keywords. That is, each 
variable in the variable list has the form numeric variable [“ dim 1[, dim 2]>], and dim 1 and dim 2 are 
real numeric expressions. 





Not usable in CALC mode. 
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COMPLE and COMPLE SHORT both allocate memory for variables and arrays. If the array or 
variable does not already exist, creation occurs upon execution of the COMPLE or 
COMPLE? HORT statement, and all variables and array elements are initialized to (0,0). The 
dimension limits of arrays are evaluated at creation time. The lowest numbered subscript in any 
dimension is 0 or 1, depending upon the OF TIOH BASE setting when the array is created. 





A COMPLES statement redimensions existing arrays if they are type COMPLE, but does not 
reinitialize them to (0,0). Similarly, a COMPLE SHORT statement redimensions existing arrays if 
they are type COMPLEX SHORT, but does not reinitialize them to (0,0). If an array is being ex- 
panded, then all newly-created elements will be initialized. Redimensioning does preserve the sequence 
of elements within an array, but not necessarily the elements’ positions within an array. Refer to the 














HP-71 Owner's Manual, section 3, under the topic Declaring Arrays (GIM, REAL, ORT, 
IHTEGER), for more information. 
The following table indicates the conditions that apply to COMPLE and COMPL SHORT vari- 





ables and arrays. 


COMPLES and COMPLE SHORT Numeric Variables 


Initial value (0, 0) 
Numeric precision 
COMPLES 12 decimal digits 
COMPLEX SHORT 5 decimal digits 
Exponent range +499 
Maximum number of array dimensions 2 
Maximum dimension limit 65535 


Simple variable memory usage (bytes) 





COMPLEX 25.5 
COMPLEX SHORT 18.5 
Array memory usage (bytes) 
COMPLE 16 = (dim 1 — option base + 1) 
+ (dim 2 — option base + 1) + 9.5 
COMPLEX SHORT 9 + (dim 1 — option base + 1) 


+ (dim 2 — option base + 1) + 9.5 
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Complex Number Operations 


(,) Real to Complex Conversion 


<X,Y) 
where X and Y are real- or complex-valued numeric expressions. 
This is the way the HP-71 recognizes a complex number: as an ordered pair of real numbers. Since (X,Y) 


is defined as (real part of X, real part of Y), if either X or Y is complex, (X,Y) is not necessarily equivalent 
toX + iY. 


Can be used in CALC mode. 








REPT Real Part of Complex Number 





REPT¢Z> 


where Z is a real- or complex-valued numeric expression. 





Returns the real part (first component) of Z. If Z is real, REP T*Z) = Z. 
Can be used in CALC mode. 





IMPT Imaginary Part of Complex Number 








IMNPT&Z3 


where Z is a real- or complex-valued numeric expression. 





Returns the imaginary part (second component) of Z. If Z is real, IMP T*<Z> =0. 








Can be used in CALC mode. 
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Other Complex Operations 


The Math Pac allows extension of many operations of the HP-71 to the complex case. These include 
numeric functions such as SIN, #, etc., as described in section 5. Other extensions are the ability to 
assign values to complex variables created by a COMPLEX or COMPLEX SHORT statement, execution 
of the RES function when the last result is complex, and so on. In other words, when the Math Pac 
module is plugged in, the HP-71 can operate with complex numbers in much the same way that it 
operates with real numbers. 





An important feature provided by the Math Pac is the extension of IMAGE format strings to include 
complex field specifiers. This extension is described below. Refer to the IMAGE keyword entry in the 
HP-71 Reference Manual for additional information on format strings. 


C(,) Complex Field in an IMAGE String 


| [n]o «format string» 


where rn is an optional multiplier. 


Causes a complex expression in a [:1SF or PF IHT output list to be formatted according to the format 
string. The real part is formatted first and the imaginary part second. On output, the number is enclosed 
in parentheses, with the real and imaginary parts separated by a comma. The comma is sent out when 
the second numeric field is encountered. 


| The format string may not include: 


e A Carriage control symbol (#). 
° String fields. 
« Imbedded complex format strings. 


The format string must include two and only two numeric specifiers, but no special restrictions (other 
than those stated above) are placed on non-numeric specifiers. 





Not usable in CALC mode. 
Complex expressions ina (ISP USING or FEINT USING output list may only be formatted by a 
complex field in the IMAGE list. Likewise, real expressions ina [ISP USING or PFEINT USTHG 
output list may not be formatted by a complex field in the IM 
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COMPLEX, COMPLEX SHORT, (,), REPT, IMPT 


Input/Result 


DESTROY ALL (END LINE 





SOHPLEX Z, Wit3d, VC?; 


END LINE 


COMPLE® SHORT Ct4,73, 


END LINE 





Z2 = ¢1,9aR¢2 


© [END LINE 





2) LEND LINE 








VC6,32=3 LEND LINE 














Ved, 193068, 5 > CENDLINE 
































Insures that none of the variables and arrays in 
the following statements exist. If one did exist, it 
would not be initialized to (0,0) when the vari- 
able or array declaration statement is executed. 


Creates a complex variable, a complex vector, 
and a complex matrix. The variable = and all ele- 
ments of the arrays 41 and ‘ are initialized to 
(0,0). 

Creates a complex short array and a complex 
short variable. ‘ and all elements of { are initial- 
ized to (0,0). 


Assigns the complex number 1 + 5i to 2. 


The HP-71 representation of the complex num- 
ber 1 + 5i. 


Assigns the real number = to the complex array 
element “6,5. 


Displays two array element values. 





Complex element ‘'{ 1, 1 was assigned 

©, > at its creation. Since the real number = 
was assigned to a complex element, it becomes 
the complex number «2, >. 








Assigns (1,3) to ‘’, since (1,3) is 
CREPTC1,23,REPTC3,49%. 


Displays the complex number ‘’. 


FES (END LINE Displays the value of the most recently executed 
or displayed numeric expresion, which in this case 
is complex. 











REPTCY) ; IMPT¢Y9 


Complex IMAGE Fields 


Input/Result 


5 STD @ COMPLEX Y 

10 Y=(69.14,—12.7) 

20 DISP USING 100; Y 

30 DISP USING 200; Y,Y 

40 DISP USING 300; Y,Y 

50 DISP USING 400; Y,Y,Y 

60 DISP USING “C(DDD,DDD)’;Y 
100 IMAGE C(2D.2D,4D.2D‘’’) 

200 IMAGE C(4Z,XXX,4%),/,C(4Z,XXX4*) 
300 IMAGE C(B,K'i"),X,C(*,4*.2DE) 
400 IMAGE 3C(2(DDD,xx)) 


RUN 


Line 100 IMAGE display. 
Line 200 IMAGE display. 
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Line 300 IMAGE display. 


Line 400 IMAGE display. 


Line 60 display. 





Section 4 


Real Scalar Functions 





Hyperbolic Functions 


The functions S1HH, OSH, and TAHH (described below) are also defined for complex arguments. See 
section 5. 


SINH Hyperbolic Sine 








SIWH¢X? 


where X is a real-valued numeric expression, |X| < 1151.98569368 





Can be used in CALC mode. 





COSH Hyperbolic Cosine 





COSHEX2 







where X is a real-valued numeric expression, |x| < 1151.98569368 





Can be used in CALC mode. 








TANH Hyperbolic Tangent 





TAHHEX> 


where X is a real-valued numeric expression. 





Can be used in CALC mode. 
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ASINH Inverse Hyperbolic Sine 





ASIHHEX? 


where X is a real-valued numeric expression. 





Can be used in CALC mode. 





ACOSH Inverse Hyperbolic Cosine 





| ACOSHEX? 


where X is a real-valued numeric expression, X > 1. 





Can be used in CALC mode. 





} 
| 
| 
l 





ATANH Inverse Hyperbolic Tangent 





ATAHH«X> 


where X is a real-valued numeric expression, —1 < X < 1. 








Can be used in CALC mode. 








Other Functions Performing Calculations 


GAMMA Gamma Function 





| 
GAMMALXS 


where X is a real-valued numeric expression whose range is defined as follows: 
X not equal to zero or a negative integer. 
—253 < X < 254.1190554375. 


Within the range —263 < X < —253, certain values of X cause GAMNMHA*X2> to underflow as indicated 
by the graph of GAMMAEX>. 





For X < —263, |GAMMA“X3| <.MINREAL, so GAMMAX> will always underflow here. 
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GAMMA (continued) 





If X equals a positive integer, GANMACX > = FACTEX—13. 
In general, GAMMA®X2 = T (X), defined for X > Oas 


re = J ee-tat 


and defined for other values of X by analytic continuation. 


Can be used in CALC mode. 





LOG2 Base 2 Logarithm 


Lo 





where X is a real-valued numeric expression, X > 0. 





LOG2¢X) = log,(x) = ae 


Can be used in CALC mode. 





SCALE10 Power of Ten Scaling 





SCALE1L 


AEX, P? 


where X is a real-valued numeric expression and P is a real numeric expression that must evaluate to an 


integer value. 
k 


Multiplies X by 10 raised to the power P by adding P to the exponent of X. You will find SCALEI4 
useful in preventing intermediate underflows and overflows in long chain calculations. 


Can be used in CALC mode. 
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Integer Round 


IROUND Round to Integer 





TROUHOGX? 


where X is a real-valued numeric expression. 





Rounds X to an integer using the current OP TION ROUND setting. | 


| Can be used in CALC mode. 


Functions Providing Information 


NANS$ Not-a-Number Diagnostic Information 





HANS CX? 


where X is a real-valued numeric expression. 





Returns a string representing the error number contained in its !4s' argument; that is, the number of 
the error that caused the =! to be created. The string returned is of the same form as the number 
returned by the ER RM function (refer to the HP-71 Reference Manual). However, the LEX identification 
number is 0 for all a's created by Math Pac functions since the Math Pac uses only HP-71 error 
messages when creating ls. 


If X is not a Hatt, then HAN#«X> returns a null string. 


Not usable in CALC mode. 





NEIGHBOR Nearest Machine Number 





HEIGHBOREX,Y3 


where X and Y are real-valued numeric expressions. 
Returns the nearest machine-representable number to X in the direction toward Y. This is the machine 
successor (or predecessor) of X depending on the relative location of Y. You will find HE 1GHEOF useful 
when you wish to evaluate a function in a local neighborhood of a given value. 





| Can be used in CALC mode. 











TYPE 
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Expression Type and Dimension 





TYPE*X) 


where X is a real-, complex-, string-, or array-valued 


expression. 





Returns an integer from 0 through 8 depending on the type and dimension of X as shown in the follow- 


ing table. 


Except for string and array arguments, can be used 





in CALC mode. 








x 


TYPE CX) 





variables.) 


variables.) 

Simple string 

INTEGER array 

SHORT array 

REAL array 

COMPLEX SHORT array 
COMPLEX array 





String array 





Simple real (includes INTEGER, 
SHORT, and REAL simple 


Simple complex (includes COMPLE 
and COMPLEX SHORT simple 


cm 


i) fae 


fe 


an 


Now om 











Examples 


COSH, SINH, ATANH, ACOSH 
Input/Result 


3 [END LINE 














Hyperbolic cosine of a numeric constant. 
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SINHCis2te2 








END LiNE Hyperbolic sine of a numeric expression. 














END LINE Inverse hyperbolic tangent of a numeric ex- 
ATAHHe 1. ; 3 pression with a numeric variable. 














>> [END LINE Inverse hyperbolic cosine of a numeric 
expression. 





LOG2, IROUND 





Input/Result 
(END LINE } Logarithm (base 2) of a numeric expression. 
OPTION ROUND NEAR (ENDLINE) Rounds to the nearest integer (the nearest even 


integer in case of a tie). 


TROUHD ¢ 





34, S > [END LINE 











OPTION ROUND POS (ENDLINE Rounds to the nearest larger integer. 


TROUMO’ 234,59 (ENDLINE 
































Section 4: Real Scalar Functions 33 





NANS, NEIGHBOR, TYPE 











Input/Result 
H=TRAPCIVL, 2> (ENDLINE Sets trap value = for I\'L. Refer to the HP-71 
Reference Manual for information on the TEAF 
function. 
IH «IHF > [END LINE 
WRH: Inwalid Ara Trap value 2 for I‘ causes a warning, not an 





error, to be given when the invalid operation 
SIM« LHF» is executed. 





* [END LINE 











The invalid operation assigns 'at/ (Not-a- 
Number) to *, since I'L has a trap value of =. 











END LINE 


il The message number associated with the value 
Hal identifies the Invalid Ara message. 


| HEIGHBOR« 1,5 > (END LINE 








The nearest machine number to 1 in the direc- 
tion toward 5. 





NEIGHBOR¢1,-16 


m 
Z| 
5} 
ey 
z 
Mm, 





The nearest machine number to | in the direc- 
tion toward -18, 











NHEIGHBOR¢ 1LE468, 1E461> [ENDLINE 


The nearest machine number to 1 £48 in the 
direction toward 1E481. 
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HEIGHBORS 1 


Ai) (END LINE 








The nearest machine number to |. : 
in the direction toward 4&. | 








INTEG 3,9) 
COMPLE u 
TYPEC23; TYPECI3; TYPEC TYPE 








3 TYPECW> (END LINE 











The numbers returned by T''F'E identify the | 
type and dimension of each of the expressions. 

















Section 5 





Complex Functions and Operations 


Many useful functions are defined for complex as well as real arguments. The Math Pac allows you to 
use many HP-71 keywords for both complex and real arguments. In addition, this section describes 
other keywords defined specifically for complex operations. 


All the functions and operations described in this section (except RES, AF.G, CON, and the relational 
operators) return a complex-type result. 


With the exception of the FECT function, all complex numbers Z and W are assumed to be in rectan- 
gular, not polar, form. 


The two-dimensional nature of these functions precludes giving simple bounds for the arguments that 
will avoid underflow and overflow messages. 


Operators 


om Addition 





Z+W 


where Z and/or W are complex-valued numeric expressions. 





Can be used in CALC mode. 





al Unary Minus 





mz 


where Z is a complex-valued numeric expression. 








Can be used in CALC mode. 
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Subtraction 





Z-W 
where Z and/or W are complex-valued numeric expressions. 


Can be used in CALC mode. 




















* Multiplication 
ZW 
where Z and/or W are complex-valued numeric expressions. 
Can be used in CALC mode. 

i Division 





ZW 


where Z and/or W are complex-valued numeric expressions, W # (0,0). 


Can be used in CALC mode. 











Exponentiation 





Zw 


where Z and/or W are complex-valued numeric expressions. 


Returns the principal value of Z” = eI) 
Can be used in CALC mode. 
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Logarithmic Functions 


LOG 


Natural Logarithm 










LOGeZ> or LH&Z3 

where Z is a complex-valued numeric expression, Z # (0,0). 

If Z = x + iy, and R (cos @ + i sin 4) is the polar representation of Z, then 
LOG<Z3 =InR + i6. 

where —x < @ < 7m (radian measure). 

Can be used in CALC mode. 





EXP 


EXP<Z) 







where Z is a complex-valued numeric expression. 
IfZ =x + iy, then 
EXP¢Z) = e+ = e* (cosy + isin y). 
where y is taken to be radian measure. 
Can be used in CALC mode. 


Exponential 





All trigonometric calculations take their arguments to be in radian measure regardless of the angular 


setting. 


SIN 
SIW¢Z> 
where Z is a complex-valued numeric expression. 


lf Z = x + iy, then 


SIHeZ> = sin (x + iy) = sin x cosh y + i cos x sinh y. 


Can be used in CALC mode. 


where Z is a complex-valued numeric expression. 


lf Z = x + iy, then 


COoS¢Z3 = cos (x + ly) = cos x cosh y — j sin x sinh y. 


Can be used in CALC mode. 


TAN 
TAHEZ2 
where Z is a complex-valued numeric expression. 


If Z = x + iy, then 


TANSZ2 = tan (x + iy) 


Can be used in CALC mode. 


sin (x + jy) _ sin x cos x + i sinh y cosh y 
cos (x + iy) sinh*y + cos°x 


Sine 


Cosine 


Tangent 
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Hyperbolic Sine 





SINH¢Z2 
where Z is a complex-valued numeric expression. 


IfZ =x + iy, then 
SIHH*Z3 = sinh(x + iy) = (—/) sin (—y + ix). 
Can be used in CALC mode. 


Hyperbolic Cosine 
COSH*Z> 
where Z is a complex-valued numeric expression. 


If Z = x + iy, then 
COSHEZs = cosh (x + iy) = cos (—y + ix). 
Can be used in CALC mode. 





TANH Hyperbolic Tangent 


TAMHEZ> 


where Z is a complex-valued numeric expression. 








If Z = x + iy, then 
TANH¢Z> = tanh (x + iy) = (—/) tan (—y + ix). 
Can be used in CALC mode. 








Polar/Rectangular Conversions 


POLAR Rectangular to Polar Conversion 
POLAR CZ> 
where Z is a real- or complex-valued numeric expression. 


If Z = x + iy, and R (cos @ + i sin 4) is the polar representation of Z, then 
POLARCZ> = (R, 0) 


The angle @ is expressed in degrees (—180 < @ < 180) or radians (—x < @ < 7) according to the 
current angular setting. 


Can be used in CALC mode. 


RECT Polar to Rectangular Conversion 
RECTéZy 
where Z is a real- or complex-valued numeric expression. 


RECT is the only keyword in this section that assumes its argument Z to be in polar form. 

If Z = (R,6), where R (cos @ + / sin @) is the polar representation of the complex number x + iy, then 
RECT¢Z3 =x + ly 

The angle 4 is taken to be in degrees or radians according to the current angular setting. 

Can be used in CALC mode. 


SQRT Square Root 
SHRTCZ2 or SUR¢Z3 
where Z is a complex-valued numeric expression. 


Returns the complex principal value of the square root of Z. 


Can be used in CALC mode. 
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Unit Vector 

| SGH¢Z> 

where Z is a complex-valued numeric expression. 
_ Returns the unit vector in the direction of Z; that is, 

; 7 z x + iy 
SGH{Z) = ———_= 
|x + iy| Vx? + y? 

where Z = x + iy. 

If Z = (0,0), then SGH&Z> = Z. 

Can be used in CALC mode. 
| ed 


Absolute Value 


| ABS«zZ> 





where Z is a complex-valued numeric expression. 
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ARG Argument 





ARGEZ3 


where Z is a real- or complex-valued numeric expression. 





If Z = x + iy and R (cos @ + i sin 8) is the polar representation of Z, then 
ARG(Z3 = 90. 


The angle 4 is expressed in degrees (—180 < # < 180) or radians (—x < # < x) according to the 
current angular setting. 


His ¢Z. always returns real type. 
Can be used in CALC mode. 


CONJ Complex Conjugate 





COMIC Z% 
where Z is a real- or complex-valued numeric expression. 
If Z = x + ly, then 
CONUCZ) =x — ily 
COMJ«Z. always returns the same type (real or complex) as Z. 


Can be used in CALC mode. 








PROJ Projective Infinity 





PROIGZ3 
where Z is a real- or complex-valued numeric expression. 
If Z = x + iy, then 
PROIJEZ>s = Z if ABS¢Z3 # Inf 


or 


>= Inf + if 





Can be used in CALC mode. 
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Relational Operators 


a, >, #, ? Equal or Unordered 
















Z comparison operator W 


where Z and/or W are complex-valued numeric expressions. 





When at least one of two expressions is complex valued, only two comparison results are possible: 
either the expressions are equal or they are unordered (or unequal, which is equivalent to unordered in 
this case). 


| Suppose Z = x + iy andW =u + Wv. 

\fx = u and y = v, then any comparison that contains = is true (that is, evaluates to 1). 
Ifx # u ory # v, then any comparison that contains # or ® is true. 

Any comparison that contains « or » without ? or # produces an exception. 


Can be used in CALC mode. 











Examples 


ame *, / 


Input/Result 


STD @ COMPLEX 
2=(4,5) @ 
2 +h) [END LINE] 





Z-W 
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Ce 4, 39 [END LINE 














> (END LINE 





*, LOG, EXP 


Input/Result 


END LINE 
oe, 39 (END LINE 











2 LEND LINE 








ExXPC¢1, 233 [ENDLINE 








SIN, TAN, COSH 


Input/Result 


FI%4 [END LINE] 
SINC¢21, 233 





TANCES, 539 











ABS, ARG, CONJ, PROJ 


Input/Result 
FIX4 
ABS 45 








EES [END LINE 





























END LINE 
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The fourth quadrant angle @ measured in radi- 
ans, which is the argument of the complex num- 
ber 3 — 7i. 
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STO @ CONJCe1, 





PROJC¢ 








END LINE 





POLAR, RECT, SGN 








Input/Result 

STD [END LINE] 

DEGREES 

POLAR¢-12 Rectangular to polar conversion for a real 
argument. 


The absolute value (r) is 1 and the argument (4) 
is 126 degrees. 





Fine 


FOLARC £3,433 Rectangular to polar conversion for a complex 
argument. 


= 


The absolute value (r) is 5. & 
ment (4) is 52. 1341 degrees. 





‘+ and the argu- 








RAOQIAHS [END LINE 











RECTC&-5,PI“43 > (ENDLINE Polar to rectangular conversion for a complex 
argument. The absolute value (r) is 5 and the 
argument (#) is —3m/4 radians. Since the R 
given is negative, this is the reflection of the polar 
point (5,PI/4) through the origin. 










0) 


ar 





SORT((-16,-03 


. (8,9900,-4,9998) 


LOG({-EXP¢S3, 803 
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The real part (x) and the imaginary part (y) are 
both -2,52355. 


48 Section 5: Complex Functions and Operations 


Additional Information 


In general, the inverse of a function f(z)—denoted f—1(z)—has more than one value for any argument z. 
However, the Math Pac calculates the single principal value, which lies in the part of the range defined 
as the principal branch of the inverse function f~ !(z). 


The illustrations that follow show the principal branches that the Math Pac uses for S&F T and LOG. 
The left-hand graph in each figure represents the cut domain of the inverse function; the right-hand 
graph shows the range of the principal branch. The blue and the black lines in the left-hand graph are 
mapped, under the inverse function, to the corresponding blue and black lines in the right-hand graph. 


SQRT 


Vz =Vr el"? for -x <0< 3 
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LN(z) = Inr + if for —7 <0 <4 





| 
MINI CLL 


w = LN(z) 
principal branch of w* is derived from that of the log function and the equation: 


w? = exp (z LN w), 
where LN denotes the single-valued function. 


‘To determine all values of the inverse function, use the expressions below to derive these values from 
the principal values calculated by the Math Pac. In these expressions, k = 0, +1, +2, and so on, and 
uppercase letters denote single-valued functions. 


yz = +SQR(z) In(z) = LN(z) + 2zik 


w? = werrike 


Section 6 


Array Input and Output 














The keywords in this section enable you to: 
« Fill an array with values. 
« Display or print values already in an array. 


Assignments 


= Simple Assignment 





MAT A=B 


where A and B are both vectors or both matrices. 

Array B may be either real or complex type. 

If B is complex, then A must be complex. 

If B is real, then A may be real or complex; if complex, all imaginary parts of all elements in A are set to 


zero. 















implicitly redimensions A to be the same size as B, and assigns the value of every element in B to the 
corresponding element in A. 


To halt operation, press twice. 
Not usable in CALC mode. 
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= () Numeric Expression Assignment 





MAT A=tX) 


where X is either a real- or complex-valued numeric expression. 
| If X is complex, then array A must be complex type. 
If X is real, then A may be real or complex; if complex, all imaginary parts of all elements in A are set to 


zero. 


Assigns X to all elements of A. Array A is not redimensioned. 


| To halt operation, press twice. 
Not usable in CALC mode. 


CON Constant Array 


MAT A=CON [¢X[.Y] 2] 


where A is either a real- or complex-type array, and the optional redimensioning subscripts X and Y are 
real-valued numeric expressions. X and Y are rounded to the nearest integer just as are subscripts in 
CIM statements. 


Assigns the real value one to all elements of A. If redimensioning subscript(s) are provided, A is explicitly 
redimensioned according to the number and value of those subcripts. 








Not usable in CALC mode. 





IDN Identity Matrix 


| 
| MAT A=IOW [¢X, Yo] 


where A is a real- or complex-type array and where the optional redimensioning subscripts X and Y are 
real-valued numeric expressions with the same rounded integer value. X and Y are rounded to the 


nearest integer just as are subscripts in 1! statements. If X and Y are not provided, A must be a 
square matrix (it must have two equal subscripts). 





!f no redimensioning subscripts X and Y are provided, then A will become an identity matrix. If 
redimensioning subscripts X and Y are provided, then A is explicitly redimensioned to a square matrix 
with the upper bound of each subscript equal to the rounded integer value of X and Y and then assigned 
the values of an identity matrix. 


Not usable in CALC mode. 
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Zero Array 










MAT A=ZER [¢X[.¥] 2] or NAT A=ZERO [¢X [.Y] 3] 


where A is either a real- or complex-type array, and the optional redimensioning subscripts X and Y are 
real-valued numeric expressions. X and Y are rounded to the nearest integer just as are subscripts in 
OTM statements. 





Assigns zero to all elements of A. If redimensioning subscript(s) are provided, A is explicitly 
redimensioned according to the number and value of those subcripts. 


Not usable in CALC mode. 








Array Input 


INPUT Assign Values from Keyboard Input 





MAT INPUT A[, Bl... 


where A (and B) are real- or complex-type array(s). 





Assigns real or complex numbers to the specified array(s). Complex values cannot be assigned to real 
array elements. MAT INPUT prompts with the name of an array element and then accepts a numeric 
expression from the keyboard, evaluates that expression, and assigns the result as the value of that 
element. For each array, AT INPUT gives prompts for the elements in row order (from left to right in 
each row, from the first row to the last). If there is more than one array, they are handled in the order 
specified. 

When the name of an array element is displayed, enter its value by typing in the numeric expression and 
then pressing [END LINE]. You can enter values for several consecutive elements by separating the val- 
ues with commas. When an array is filled, the remaining values are automatically entered into the next 
array. After you press [END LINE], the computer will display the name of the next element (if any) to be 
assigned a value. 
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INPUT (continued) 


In other respects, MAT IHFLT acts as does INFLUT. For instance: 


* The Command Stack is always active during MAT IHFLT execution. You can move up and down 
in the Command Stack with (A), [¥), (9) (4), and (8) Cy) without first pressing [S) [CDs]. 

* You can use a direct execute user-defined key to provide the response to the MAT IHPLIT 
prompt. 


* The [f] (VIEW) key sequence and the [9] key sequence are active during MAT IPL 
execution. 





























* If you are making a response to a NAT INFLT statement, but have not pressed [END LINE], 
pressing once clears the typed entry, allowing another entry to be typed. If you press [ATTN 
twice, the HP-71 clears the entry, pauses the program, and clears the display. 











Not usable in CALC mode. 


Array Output 


To halt the operation of any of the keywords described below you need press only once. 


DISP Display in Standard Format 


eel | C 


where A (and B) are real- or complex-type array(s). 


Displays the values of the elements of the specified arrays. The values are displayed in row order. Each 
row begins on a new line; a blank line is displayed between the last row of an array and the first row of | 
the next array. 


The choice of terminator—comma or semicolon—determines the spacing between the elements of an 
array. 


Terminator Spacing Between Elements 


Close: Elements are separated by two spaces. A minus sign, if present, 
occupies one of the two spaces. 
Wide: Elements are placed in 21-column fields. 


If the last array specified doesn’t have a terminator, the array will be displayed with wide spacing be- 
tween elements. 





Not usable in CALC mode. 
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PRINT Print in Standard Format 


Bon ELL 


where A (and B) are real- or complex-type array(s). 








Prints the values of the specified arrays. Operation is identical to MAT [ISF, except that the output is 
sent to the PFRIHTER I5 device, which requires HP-IL. If no PR INTER 15 device is present, out- 
put is sent to the display, or to the HP-IL DISPLAY IS device. Also, you can override the CR/LF 
normally generated by MAT FRINT with the ENOLINE statement. EHODLINE is described in the 
HP-71 Reference Manual and in section 13 of the HP-71 Owner's Manual. 


Not usable in CALC mode. 





DISP USING Display Using Custom Format 


format string ’ , 
MAT DISP USING iA B |... 
line number ; ; 


where A (and B) are real- or complex-type array(s). 








Displays the values of the elements of the specified arrays in a format determined by the format string or 
by the specified 1 MAGE statement identified by the /ine number. (Refer to the HP-71 Reference Manual 
for information about DISF USIHG, format strings, IMAGE statements, and their results). 


If any array is complex type, the corresponding field specifier in the format string or IMAGE statement 
must be a complex field specifier. Refer to the description of the complex field specifier (2°. >) in 
section 3, page 22. 


The values are displayed in row order. Each row begins on a new line; a blank line is displayed between 
the last row of an array and the first row of the next array. 


The terminators between the arrays—commas or semicolons—serve only to separate the arrays and 
have no effect on the display format. 


The Math Pac must be plugged in to REHMLINEER a program containing a MAT OSPF USING [ine 
number] statement; otherwise, the /ine number will not be correctly updated. 


Not usable in CALC mode. 
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PRINT USING Print Using Custom Format 


format string Z Py 
MAT PRINT USING iA B fics 
line number i : 


where A (and B) are real- or complex-type array(s). 








Operation is identical to MAT OISP USING, except that the output is sent to the PRINTER [5 
device which requires HP-IL. If no FR IMTER 15 device is present, output is sent to the display, or to 
the HP-IL DISPLAY IS device. Also, you can override the CR/LF normally generated by 
MAT PRIHT USING with the ENDLINE statement. ENDLINE is described in the HP-71 Ref- 
erence Manual and in section 13 of the HP-71 Owner's Manual. 


Not usable in CALC mode. 














Examples 


With the optional delay of & or larger (infinite line replacement delay), you press (or any 
other key) to display the next line. So you can control how long each array row is displayed. 


CON, IDN, ZER, DISP 


Input/Result 

OPTION BASE 1 @ STD 

DIM ACS, 3), BC 1> B is dimensioned to be a one element vector. 
COMPLEX C¢18, 26% 

MAT ASTON 


MAT. (ESP Ay Displays the identity matrix A with close spacing 
between the elements. 














MAT B=ZERC2 





» 23 (END LINE 


MAT ISP B: 





MAT C=CONC3, 3) 


MAT DISF (; (ENDLINE 








INPUT 


Input/Result 

OPTION BASE 1 (ENDLINE) 
DIM A¢2,33,B¢3> CENDUNE) 
OPTION BASE & (ENDLINE) 
COMPLEX C¢2, 13 (ENDUINE) 









NAT INPUT A, E, C [ENDLINE) 
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Redimensions B from a one-element vector to a 
2 x 2 matrix and assigns to it a zero array. 


Redimensions C and assigns to it a constant 
array. 


Declares C to be a 3 x 2 complex array (remem- 
ber we are in OF TION BASE &), 


Prompts for the first element’s value. 


More than one value can be entered. 


Prompts for the fifth element’s value. 
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5,6, 7 [ENDLINE 
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» 1 [END LINE] 





HAN (END LINE) 
STO @ MAT OISF 


I he 


2 








2 (END LINE 


A; 


B 








Enters values for the last two 
elements of A and the first element of B. 


Enters values for the last two elements of B and 
the first element of the complex array C. 


Enters values for the next four elements of C. 


Enters “not a number” for the last element of C. 


Displays each array in sequence, with a blank 
line between each. 


and 


f C. 


DISP USING 


Input/Result 


10 OPTION BASE 1 @ INTEGER A(5,5) 
15 WIDTH 22 @ DELAY 8 


20 COMPLEX SHORT Z(3,4) 
25 MAT A=IDN @ MAT Z=((4,5)) 


30 MAT DISP USING ‘DDD,ZZZ';A,A 


35 MAT DISP USING ‘#,D';A @ DISP 4 


40 MAT DISP USING 100;Z 
45 DELAY 1 
100 IMAGE C(K,2D,'i’) 
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Causes the output to appear in the display as 
shown below. After each display, press [END LINE 
to produce next display. 





Assigns the identity matrix to A and the 


c 


complex number ‘4, 5% to every element of Z. 
This format string consists of two field 


specifiers, DUD and 222. Each element of A is 
displayed according to these field specifiers used 
repeatedly until all elements have been displayed. 
The final element of A is displayed according to 
O00. Then a blank line is displayed, followed by 
another display of all elements of A. The field 
specifier 222 (the next specifier in the format 
string) is used to format the display of the first 
element during this second display of A. 


The # symbol supresses the automatic end-of- 
line sequence (CR/LF) following the display of A. 
This causes 4+ to be displayed on the same line 
as the last element of A. 





The IMAGE statement must use the C¢, » form 
to format the display of a complex array. The 
parentheses must contain two numeric field 
specifiers. 
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RUN 


The O format symbol replaces leading zeros with 
blanks. Since A is an identity matrix, element 
(1,1) is 1. Therefore the two leading zeros are 
replaced with blanks, and element (1,1) is dis- 
played as i. The 2 format symbol fills each 
leading zero with 4, so element (1,2) is displayed 
as 444. The remaining elements, in row order, 
d according to the format string 


= used repeatedly. 








After the last (fifth) element of the first row is 
displayed, an end-of-line sequence (carriage re- 
turn, line feed) is sent, causing the display of 
element (2,1) to start a new line. 


The field specifier DOC formats the display of 
the last element of A, causing the display of 1. 


Following the display of the last element of the 
last row, a second end-of-line sequence is sent, 

causing the display of a blank line between the 
two displays of array A. 


Since the variable list following the format 
string in line 30 is A, A, array A is displayed 
twice. This time, element (1,1) is displayed 
according to the field specifier , since OOO 
was used just above for the last element of A 
during the first display of this array. 








Since this is the display of the last array in the 
variable list of line 30, no blank line is displayed, 
even though this display line ends with the last 
element of the last row of A. 
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Since the portion of the format string of line 35 
that controls character display consists only of 0, 
the elements of each row of A are displayed with 
no extra characters or spaces. 


The # symbol in the format string of line 35 
supresses the end-of-line sequence normally sent 
after the display of the final row of the last array 
in the variable list. 


The symbol F: in the format string of line 100 
specifies a compact field, resulting in the display 
of no leading or trailing blanks. This symbol 
controls the display format of the real part of 
each (identical) element of Z. The display of the 
imaginary part of each element is controlled by 
20), Since the imaginary part, 3, consists of only 
one digit, a leading blank is displayed. The com- 
plex image specification [* * causes the display 
of the parentheses and comma. 





The display of each row is ended with an end-of- 
line sequence, so each new row starts a new dis- 
play line. 


Section 7 


Array Arithmetic 





The keywords in this section perform arithmetic operations on arrays. The dimensions of the operand 
arrays must be compatible with the particular operation, as discussed below. 


For addition and subtraction, the operand arrays must both be vectors or both be matrices, and 
they must have the same number of rows and the same number of columns. In this case we will say 
that the arrays are conformable for addition. 

* For multiplication of two arrays, the first array must be a matrix, while the second array can be a ° 
matrix or a vector. The number of columns of the first array must be equal to the number of rows 
of the second array. If these conditions are satisfied, we will say that the arrays are conformable for 
multiplication. 

* For transpose multiplication of two arrays, the first array must be a matrix, while the second array 
can be a matrix or a vector. The number of rows of the first array must be equal to the number of 
rows of the second array. If these conditions are satisfied, we will say that the arrays are con- 
formable for transpose multiplication. 





Operators 
=— Negation 
MAT A=-B 


where A and B are both vectors or both matrices. 

Array B may be either real or complex type. 

If B is complex, then A must be complex. 

If B is real, then A may be real or complex; if complex, all imaginary parts of all elements in A are set to 
zero. 





Implicitly redimensions A to be the same size as B and assigns to each element of A the negative of the 
corresponding element of B. 


To halt operation, press [ATTN] twice. 








Not usable in CALC mode. 
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oe Addition 





NAT A=B+C 


where A, B, and C are all vectors or all matrices, and B and C are conformable for addition. 

Arrays B and C may be either real or complex type. 

lf either B or C is complex, then A must be complex. 

If both B and C are real, then A may be real or complex; if complex, all imaginary parts of all elements in 
A are set to zero. 





implicitly redimensions A to be the same size as B and C, and assigns to each element of A the sum of 
the values of the corresponding elements of B and C. 


To halt operation, press [ATTN] twice. 





Not usable in CALC mode. 





= Subtraction 





MAT A=B-C 


where A, B, and C are all vectors or all matrices, and B and C are conformable for addition. 

Arrays B and C may be either real or complex type. 

If either B or C is complex, then A must be complex. 

If both B and C are real, then A may be real or complex; if complex, all imaginary parts of all elements in 
A are set to zero. 





implicitly redimensions A to be the same size as B and C, and assigns to each element of A the dif- 
ference of the values of the corresponding elements of B and C. 


To halt operation, press [ATTN] twice. 








Not usable in CALC mode. 











—_ 
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( )* Multiplication by a Scalar 





MAT A=*X3 +B 


where A and B are both vectors or both matrices and X is a numeric expression. 

Array B may be either real or complex type and expression X may be either real or complex valued. 

If either B or X is complex, then A must be complex. 

If both B and X are real, then A may be real or complex; if complex, all imaginary parts of all elements in 
A are set to zero. 





implicitly redimensions A to be the same size as B and assigns to each element of A the product of the 
value of X and the value of the corresponding element of B. 


To halt operation, press twice. 
Not usable in CALC mode. 





* Matrix Multiplication 


MAT A=BtC 


where B is a matrix, A and C are both vectors or both matrices, and B and C are conformable for 
multiplication. 

Arrays B and C may be either real or complex type. 

If either B or C is complex, then A must be complex. 

If both B and C are real, then A may be real or complex; if complex, all imaginary parts of all elements in 
A are set to zero. 





Implicitly redimensions A to have the same number of rows as B and the same number of columns as C. 
The values of the elements of A are determined by the usual rules of matrix multiplication. 


To halt operation, press [ATTN | twice. 





Not usable in CALC mode. 
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TRN xX Transpose Multiplication 


| MAT A= TRH CBS#C 

| where B is a matrix, A and C are both vectors or both matrices, and B and C are conformable for 
transpose multiplication. 
Arrays B and C may be either real or complex type. 
If either B or C is complex, then A must be complex. 
If both B and C are real, then A may be real or complex; if complex, all imaginary parts of all elements in 
A are set to zero. 





implicitly redimensions A to have the same number of rows as the number of columns in B and the same 
| number of columns as C. 


The result of this operation is the same as if the transpose of B (or the conjugate transpose of B, if Bis 
complex type) was computed and then post-multiplied by C. However, the Math Pac uses special mul- 
tiplication rules so that B does not have to be explicitly transposed prior to the multiplication. 


To halt operation, press [ATTN] twice. 








Not usable in CALC mode. 





Examples 


+, *, ()*, TRN * 

Input/Result 

OPTION BASE 1 @ STO (ENDLINE] 
? 43 

















% 9,E¢9) 
MAT A=ION¢E 
MAT C=¢¢3,4934A C is redimensioned to 2 x 2 and every element 


of C is assigned the product of the complex num- 
ber “3,42 and the corresponding 
element of A. 


MAT OISF C; LENDLINE 





The array C. 





or | 


MAT A=COH @ MAT C=C+A (ENDLINE 





MAT OISP C; [ENDLINE 








MAT B=A+A 


MAT DISF B; 





MAT INPUT o 





1,29,¢3,4) 
MAT E=TRA¢C > +0 [ENDLINE] 


MAT OISP € [ENDLINE] 
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C holds the array sum of A and C. No 
redimensioning is necessary since C is already the 
correct size. 


The array C. 


B is redimensioned to 2 x 2 to hold the matrix 
product A#A. 


The array B. 


E is redimensioned to be a 2 element vector to 
hold the product of the conjugate transpose of C 
and the vector D. 


The array E. 





Section 8 


Scalar-Valued Array Functions 


The keywords in this section are functions that use real- or complex-type arrays as arguments (except 
GET uses only real arrays) and give a real number as a value (except [iT can give either a real or 
complex number). Like other HP-71 functions, they may be used alone or in combination with other 
functions to produce numeric expressions. 


Determinant Functions 


DET Determinant 


DETC AS 


where A is a square real-type matrix. 


Returns the determinant of the matrix A. 
To halt operation, press twice. 
Not usable in CALC mode. 


DETL Determinant of Last Matrix 


DETL or DET 


Returns the determinant of the last real-type matrix that was: 


© Inverted ina MAT... IH! statement (described in section 9). 
5 statement (described in section 9). 






e Used as the first argument of a MAT. . 


QETL retains its value (even if the HP-71 is turned off) until another MAT... IN\ (with a real type 
argument) or a MAT. . 'S (with a real type first argument) is executed. 





Not usable in CALC mode. 
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Array Norms 


CNORM 
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One-Norm (Column Norm) 





CHORMCAS 


where A is a real- or complex-type array. 


absolute value of a complex number. 





To halt operation, press [ATTN | twice. 


Not usable in CALC mode. 


Returns the maximum value (over all columns of A) of the sums of the absolute values of all elements in 
a column. Refer to the keyword description for HES, page 41 in section 5, for the definition of the 





RNORM 


Infinity Norm (Row Norm) 





RMORM CAS 


value of a complex number. 





jo halt operation, press [ATTN] twice. 


Not usable in CALC mode. 





FNORM 


where A is a real- or complex-type array. 


Returns the maximum value (over all rows of A) of the sums of the absolute values of all elements in a 
row. Refer to the keyword description for RES, page 41 in section 5, for the definition of the absolute 


Frobenius (Euclidean) Norm 





| FHORM CA? 


L = = 


complex number. 











To halt operation, press [ATTN] twice. 





Not usable in CALC mode. 


where A is a real- or complex-type array. 


Returns the square root of the sum of the squares of the absolute values of all elements of A. Refer to 
the keyword description for HES, page 41 in section 5, for the definition of the absolute value of a 








m) 
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Inner Product 


DOT Inner (Dot) Product 





DOT«X.Y> 


where X and Y are real- or complex-type vectors with the same number of elements. 





Returns XsY, the inner product of X and Y. If both X and Y are real, then the result is real. If either X or Y 
is complex, then the result is complex. 


If X is a complex vector, then the complex conjugates of the elements of X are used to compute the 
inner product. 


To halt operation, press twice. 
Not usable in CALC mode. 











Subscript Bounds 


The following functions are useful in keeping track of array option base, number of dimensions, and 
size in each dimension, since these quantities may change when variables are dimensioned and 
redimensioned. 


UBND Subscript Upper Bound 





UBMOCA,N? or UBOUNDCA,N2 


where A is a real- or complex-type array and N is a numeric expression whose rounded integer value is 1 
or 2. 









Returns the upper bound of the Nth (first or second) subscript of A. If A is a vector, WEMOCA, = 
=1. 





Not usable in CALC mode. 








2d Array Functions 





LBND Subscript Lower Bound 


LEHOCA,N® or LEOUHOCA,N? | 
| 

where A is a real- or complex-type array and N is a numeric expression whose rounded integer value is 1 

or 2, 








BASE setting in effect when A was last dimensioned. If A is a vector, 


Returns the © 
LBHOCA, 23 = 






| Not usable in CALC mode. 


Se 


Examples 


DET, DOT 


Input/Result 

OPTION BASE 1 (ENDLINE) 
OIM Aci,16> [ENDUINE) 
MAT A=10N [END LINE) 


4A (END LINE Assigns —2 to each diagonal element; all other 
elements remain zero. 











DET A® (END LINE Displays the determinant of A. 























MAT ASTON 

MAT A=C294A Assigns 2 to each diagonal element; all other 
elements remain zero. 

MAT A=INVEA® (END LINE Computes the inverse of A. 

QET (END LINE Displays the determinant of the last real matrix 
inverted in a MAT. ..IN\ statement or used as 





the first argument of a MAT... statement. 
Refer to pages 77-79 in section 9 for definitions of 
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>, BC16) 























MAT A= END LINE Assigns = to each element of A. 
MAT B=COM (END LINE Assigns one to each element of B. 
DOTCA, B» Displays the inner product of A and B. 























COMPLE C18 (END LINE 

MAT C » 29> (END LINE] Assigns the complex number «1, 22 to each 
element of C. 

DOT<C,. AD Displays the inner product (a complex number) 
of C and A. 





RNORM, CNORM, FNORM, UBND, LBND 
Input/Result 

OPTION BASE 1 

oIn > (END LINE) 
4 [END LINE} Assigns 1 to each element of A. 
» (END LINE] Displays the row norm of A. 
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4 (END LINE 





COMPLEX SHORT 


MAT INPUT A 






C1, 
C11 a 


END LINE 
RHMORM CAD 














COMPLEX BC3> 
UBHOCA, 1); UBNOCA, 











UBHOCE, 12; UBNDe 


Displays the row norm of A. 


Displays the column norm of A. 


Displays the Frobenius norm of A. 


First, displays the upper bound of A’s first sub- 
script, then displays the upper bound of A’s sec- 
ond subscript. 


First, displays the upper bound of B’s first sub- 
script, then attempts to display the upper bound 
of B’s second subscript. Since B has only one 
subscript, UBHO®B,.2> returns -1. 
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LENOCA, 19 Displays the OF TION EASE setting when A 
was last dimensioned. 





Section 9 


Inverse, Transpose and System Solution 


Operations 


INV Matrix Inverse 











MAT A=IHVC Bo 


where A is a matrix and B is a square matrix. 

Array B may be either real or complex type. 

If B is complex, then A must be complex. 

If B is real, then A may be real or complex; if complex, all imaginary parts of all elements in A are set to 


zero. 
implicitly redimensions A to be the same size as B and assigns to A the value of the matrix inverse of B. 
To halt operation, press twice. 

Not usable in CALC mode. 


TRN Matrix Transpose or Matrix Conjugate Transpose 





MAT A=TRH*¢B? 


where A and B are matrices. 

Array B may be either real or complex type. 

If B is complex, then A must be complex. 

If B is real, then A may be real or complex; if complex, all imaginary parts of all elements in A are set to 


zero. 





implicitly redimensions A to be the same size as the matrix transpose of B. If B is real, assigns to A the 
value of the matrix transpose of B. If B is complex, assigns to A the values of the matrix conjugate 


transpose of B. 
To halt operation, press twice. 





Not usable in CALC mode. 
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Solving a System of Equations 


The Math Pac provides a quick and accurate way to solve a system of linear equations involving real or 
complex coefficients. The first step in using this capability is to translate the system of equations into a 
triple of arrays: the result array, the coefficient array, and the constant array. The result array corre- 
sponds to the variables in the equations; the coefficient array holds the values of the coefficients of the 
variables; the constant array holds the values of the constants in the equations. For example, if you 
wanted to solve the system of equations 


bx + 3y + 22 = 4 
ix + y+ 32 = 14 
6x + 4y + 9z = 1 


then the result array would correspond to the array 


x 


” 


the coefficient array would be 


and the constant array would be 


14 
1 
If we denote the result array by X, the coefficient array by A, and the constant array by B, then the 


system of equations can be written in matrix notation as AX=B. This is the form assumed by the 
keyword. 
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sys System Solution 


A,B? 








MAT X 


where A is a square matrix, X and B are both vectors or both matrices, and A and B are conformable for 
multiplication. Refer to the beginning of section 7, page 63, for a definition of “conformable for 
multiplication.” 


| Arrays A and B may be either real or complex type. 
If either A or B is complex, then X must be complex. 


If both A and B are real, then X may be real or complex; if complex, all imaginary parts of all elements in 
X are set to zero. 


Implicitly redimensions X to be the same size as B and assigns to X the computed solution to the matrix 
equation AX=B. 


To halt operation, press twice. 








Not usable in CALC mode 





Examples 


INV, TRN 


Input/Result 

OPTION BASE 1 (ENDLINE) 
DIM ACB, 39 
MAT A=I0N [ENDLINE) 


MAT | 204A (END LINE Assigns = to all diagonal elements of A. All other 
elements are zero. 





} 


xa 
u 





MAT A=IMVCAS 





END LINE 








MAT OISP A; [ENDLINE] Displays the inverse of A. 
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OIM c END LINE 
MAT C 
MAT De; 
1 
1 
1 
DIM O«¢ 
MAT D=TRNEED 


MAT DISP OD; 





i a 
4. A 
COMPLEN SHORT O¢2, 


MAT D=¢e1, 299 


MAT OISP Oj 





MAT D=TRNCOD 


MAT DISP O; 





Assigns one to all elements of C. 


Displays C 


Computes the transpose of C and redimensions 
D to be a 2 X 3 matrix. 


Displays the transpose of C. 


Assigns the complex value «i, 2° to all ele- 
ments of D. 


The complex matrix D. 


Redimensions D to 3 x 2 and assigns D the 
value of its conjugate transpose. 


The conjugate transpose of D. 
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MAT IHPUT © (ENDLINE 








@, 18> [END LINE 








€1,1> 
MAT O 








The complex matrix C. 





MAT D=INWEt> Redimensions D to 3 x 3 and assigns to D the 
value of the matrix inverse of C. 


MAT OISP OD; 


‘The inverse of the complex matrix C is the 
matrix 
10+i —24+6i —3-—2i 
9-31 8i —3-2i 
A 





—1-2i 





—24+2i 
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SYS 

To solve the system of equations given on page 78, namely, 
5x + 38y + 2z2 = 4 
Ix +y + 32 = 14 
6x + 4y + 92 = 1 

we could use the following keystrokes. 


Input/Result 

OPTION BASE 1 @ STO (ENDLINE) 
OM 3), B03>,A¢3, 23> (ENDLINE) 
MAT IMPUT 6, A (END LINE) 











4,14, 1 [END LINE Assigns values to the elements of B. 








4,9 Assigns values to the elements of A. 












(A, B® (END LINE 
END LINE Displays the values of the result array elements. 











Wd 
re HR 





Although in typical applications the result array X and constant array B are each one column arrays, 
=‘ does not restrict these arrays to only one column. This allows you, for example, to simultaneously 
solve any number of different systems, limited only by memory, of n equations in n unknowns, provided 
that the coefficients in each systems of equations are identical. The following example illustrates this 
use of 34/5 
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Example. Your company’s Publications Manager wants to determine the cost factors used by her two 
outside printers. She knows that each printer estimates jobs based on the number of pages and the 
number of photographs, plus a fixed setup charge. Given the three estimates from each printer shown 
below, write a program that calculates their cost per page, cost per photograph, and setup charge. 




















Number of | Number of Total Cost 
Job 
Pages Photographs Printer A | Printer B 
1 273 35 $5835.00 | $7362.50 
2 150 8 $3240.00 | $4085.00 
3 124 19 $2775.00 | $3517.50 








We need to solve the following system of equations for two sets of cost estimates. 
273x, + 35x + x3 = cost, 
150x, + 8x2 + x3 = costy 
124x; + 19x2 + x3 = costs 


These equations can be represented in matrix notation as AX = B, where: 


* Ais the coefficient matrix, having the number of pages in its first column, the number of photo- 
graphs in its second column, and the number of setup charges (one for each job) in its third col- 
umn. Each row contains this data for a different job. 


* B is the constant array. Each row contains cost estimates for one job from the two printers; each 
column contains one printer’s cost estimates for the three jobs. 


* X is the result array, having the unknown cost factors x,, x:, and x3 in its rows. x, is the cost per 
page, x» is the cost per photograph, and x3 is the setup charge. Since we are solving two systems, 
the constant array is a two-column matrix. So the result array must also be a matrix; that is, it 
should be declared with two dimensions. (Its size, if not the same size as that of the constant array 
B, will automatically be redimensioned to the size of B when the 2° 5 statement is executed). Each 
column will contain the cost factors for one printer. 








10 OPTION BASE 1 
20 DIM A(3,3),X(3,2),B(3,2) 


30 DATA 273,35,1 Specifications for job 1. 
40 DATA 150,8,1 Specifications for job 2. 
50 DATA 124,19,1 Specifications for job 3. 
60 DATA 5835,7362.5 Estimates for job 1. 

70 DATA 3240,4085 Estimates for job 2. 

80 DATA 2775,3517.5 Estimates for job 3. 

90 READ A,B 


100 MAT X=SYS(A,B) 


110 DISP USING ‘9A,3X,9A,/'; 
‘PRINTER A’,‘PRINTER B’ 


120 MAT DISP USING ‘X3D.2D,6X, 
3D.2D';X 
RUN 
PRIHTER A PRIHTER B 
Cost per page. 


Cost per photograph. 
Setup charge. 














Example. This example demonstrates the usefulness of ‘ in the solution of circuit analysis prob- 
lems. In the circuit shown below, the impedances of the components are indicated in complex form. We 
will determine the complex representation of the currents J; and Ip. 
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This system can be represented by the complex matrix equation 


10+200i 2001 |I7, 5 
200i (200—30)i} | 1,] | 0 


or AX =B 


Here is a program that solves for J; and I5. 


10 OPTION BASE 1 @ STD 

20 COMPLEX SHORT A(2,2),X(2) If either A or B is complex, X must be complex. 
30 DIM B(2) 

40 MAT INPUT A,B 

50 MAT X=SYS(A,B) 

60 MAT DISP X 





6>,¢68,-200),¢(6,-200>, Assigns values to the elements of A. 


1a, 3 nae 
6, 178) (END LINE} 





Assigns values to the elements of B. 
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Additional Information 





The Math Pac operations DET¢A?, MAT B=IHVCA?, and MAT X A. B®, where A is a real- 
type square matrix, all use the LU decomposition of A as an intermediary step. The method used to 
generate the LU decomposition of A is a compact Crout factorization with partial pivoting and ex- 
tended precision arithmetic. The LU decomposition of A can be represented by the equation PA = LU, 
where 


e Lis a lower triangular matrix—it has values of zero for all elements above the diagonal. 


e Uis an upper triangular martix—it has values of zero for all elements below the diagonal—with 
values of one for all elements on the diagonal. 


* Pis a permutation matrix representing the row interchanges in the matrix A resulting from partial 
pivoting. 


The factorization PA = LU is valid for any non-singular matrix A. However, special attention is paid to 
matrices A that are singular or “machine singular.” In this case, the LU decomposition is changed by an 
amount that is usually small in comparison with roundoff error. In the absence of underflow or over- 
flow, the resulting LU decomposition of A will be close, in norm, to the exact LU decomposition of 
another matrix A’, where A’ is close in norm to A. 


Consider the matrix shown below. 
1 3.0 


0 Oo 2 
.666666666667 2 0 


Although this matrix is very nearly singular, it can be successfully inverted using the IH‘! keyword: 





Input/Result 
OPTION BASE 1 (ENDLINE 
DIM ACS,33, END LINE 








MAT IHPUT A (LENDLINE 














END LINE 
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Pp ¥ 





=, & [ENDLINE A now represents the matrix given above. 





T B=IHM CAS (ENDLINE B is now the computed inverse of A. 
MAT B=B+A 
ii 





DISP 'B: Displays the identity matrix B, which is the 
product of the matrix A and its computed inverse. 





The keyword solves the matrix equation AX = B for X in several stages. First, the LU decompo- 
sition of A is found to give PA = LU. 





Using PA = LU, the equivalent problem is to solve LUX = PB for X. This is done by solving LY = PB 
for Y (forward substitution) and then solving UX = Y for X (backward substitution). This value for X is 
used as a first approximation to the desired solution in a process of iterative refinement, which pro- 
duces the final result. 


In many cases, the Math Pac will arrive at a correct solution even if the coefficient array is singular (so 
S'S to solve under- and over- 


that the formula X = A 'B is invalid). This feature allows you to use ! 
determined systems of equations. 





For an under-determined system (more variables than equations), the coefficient array will have fewer rows 
than columns. To find a solution using 5'¥ 5: 


e Append enough rows of zeros to the bottom of your coeffieient array to make it square. 


e Append corresponding rows of zeros to the constant array. 


You can now use these arrays with the 5‘'> keyword to find a solution to the original system. 





For an overdetermined system (more equations than variables), the coefficient array will have fewer columns 
than rows. To find a solution using =‘ 5: 





e Append enough columns of zeros on the right of your coefficient array to make it square. 


e Be sure that your result array is dimensioned to have at least as many rows as the new coefficient array 
has columns. 


e Add enough zeros on the bottom of your constant array to ensure conformability. 





You can now use these arrays with the keyword to find a solution to the original system. Only those 


elements in the result array that correspond to your original variables will be meaningful. 
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For both under- and overdeterminded systems the coefficient array is singular, so you should check the re- 
to see if they satisfy the original equation. 





sults returned by 5 





If A is a complex type square matrix, then MAT C=IH'!¢ A> and MAT X= A,B» use the same 
techniques as above, with the arrays A and B replaced by equivalent real-type partitioned forms. 





The 53'S keyword can also be used for inverting a square matrix A. MAT X=S'S¢A,B> will return 


the inverse of A if B is chosen to be the identity matrix. This technique is more accurate and generally 
faster than MAT X=IH\¢A®, but it requires more memory for its operation. (Refer to appendix B for 


information about memory requirements). 





Section 10 
Solving f(x) = 0 


Keywords 














You can use the keywords in this section to help you determine the solutions or minima of equations of 
from one to five real variables. 


Throughout most of this section, the operation of these keywords will be described for a one-variable 
function. Multi-variable functions are covered under the topic Nesting Rules. 


The keyword FHFOOT can be used from the keyboard or inside a program to find the value of x for 
which f(x) is zero or a minimum, provided the keyboard line or program contains the definition of the 
function. 


The keywords F\!ALUE and FIGUESS are provided to help you use FHE {OT and to interpret its re- 
sults. Since all three keywords are numeric-valued functions, they can be used alone or in combination 
with other functions and variables to form numeric expressions. A fourth keyword, F \| AR, represents 
the variable in the function being solved by FHFO0T. It also contains the most recent guess generated 
by an executing FHROOT. 


FNROOT Function Root 





FHROOTCA,B.F3 


where A, B, and F are real numeric expressions. 





Seeks a real root of the function F, starting with the two guesses A and B. These guesses can be equal, 
but if so, one is immediately perturbed. 


Returns the first value found that is one of the following: 
1. An exact root of the specified function. 
. An approximation to a root of the specified function, correct to 12 digits. 


2 
3. An approximation to a local minimum of the absolute value of the specified function. 
4. In a region where the specified function is constant. 

5 


. +£9.99999999999E499 if the search for a root led beyond the range of representable numbers. 
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FNROOT (continued) 


Not usable in CALC mode. Refer to page 97 for more information about FHF:00T and CALC mode. 


Refer to pages 97-99 for information about FNF:OOT nesting and about the interactions between 
FHROOT and [ATTN] and between FHROOT and user-defined functions. 














FVAR 


Function Variable 


FYVUAR 


Represents the variable x in f(x), the variable whose value FMFROOT seeks. | 
Also returns the most current guess generated by a running FHROOT. | 
Can be used in CALC mode. | 


FVALUE Function Value 


FVALUE 


Returns the value of the function F (the third argument of F HOUT) at the result generated by the most 
recently completed FHROOT. 


F\ALLIE retains its value, even if your HP-71 is turned off, until F W#U0T is again completed. 


Can be used in CALC mode. 


FGUESS Previous Estimate of Function Root 


Returns the next-to-last value tried as a solution in the most recently completed F Hf i0)T statement. 


FGUESS retains its value, even if your HP-71 is turned off, until FHROOT is again executed. 


Can be used in CALC mode. 
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checking the values of F'\!ALUE and FGUESS, you can interpret the result of FHEOOT as follows: 





e If FYALUE = 0, the result of FHROOT is an exact root and the result of FGUE 
number close to the root. 


will be a 


© If the result of FHROOT and the result of FGUESS differ only in the twelfth significant digit, and 
FVALUE and F< FGUESS> have opposite signs, these two numbers surround the exact root. 


e If the result of FHROOT and the result of FGUESS differ, but F'\ALUE and the value of the 
function at FGLIESS are equal, these results lie in a region where FHF is constant. 









solve an equation for a particular variable, use this procedure: 
. Write the equation to be solved in the form f(x) = 0. 


. Substitute the keyword F\!FiR everywhere for the variable you wish to solve for in the formula 
defining f(x). 


3, Use the defining formula for /(x) as the third argument for FHROOT. 


. Choose two initial guesses (which may be equal) and use these as the first two arguments for 
FNROOT, Even if only one initial guess is used, use it for both A and B, since FHROCMT always 
requires three arguments. 





Examples 


Solving x? = 2 (FHROOT, FVALUE, FVAR) 


The following six examples illustrate various ways FMF:O10T and F\\AF can be used to solve the equa- 
tion x* = 2. Initial guesses of 1 and 2 are used. The first and sixth examples show the solution, 


Example One. 


Input/Result 
FHROOTC1, 2,FVARCS 





2) FHROOT can be used from the keyboard as well 
as in a program. 


Example Two. 


10 DISP FNROOT(COS(0),LOG2(4), The initial guesses can be expressions. 
FVAR*2—2) 


20 DISP ‘FVALUE =';FVALUE 





Example Three. 


10 DEF FNG=FVAR*2—2 


20 DISP FNROOT(1,2,FNG) The third argument of FHR COT can be an 
expression or a reference to a user-defined 
function. 


30 DISP ‘FVALUE=’;FVALUE 


Example Four. 


10 DEF FNF(X)=X*2—2 
20 DISP FNROOT(1,2,FNF(FVAR)) F\'AR can appear in the user-defined function, 
as above, or in the third argument of FHROOT. 


30 DISP ‘FVALUE=’;FVALUE 


Example Five. 


10 DEF FNH The user-defined function can consist of one or 
several lines. 


20 FNH=FVAR*2—2 

30 END DEF 

40 DISP FNROOT(1,2,FNH) 
50 DISP ‘FVALUE=’';FVALUE 


Example Six. 


10 DEF FNJ(X) 
20 FNJ=X*2—2 

30 END DEF 

40 DEF FNF(X)=2*X 


50 DISP FNROOT(1,FNF(1),FNJ(FVAR)) The first or second arguments of FMR OO Tcan 
be references to user-defined functions. 


60 DISP ‘FVALUE =';FVALUE 


Input/Result 








RUN 











1.414212356235 The solution for x2 = 2. 
FVALUE = .& 
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Solving log (x) = e/x (FHEGOT, FYALUE, FVAR, FGUESS) 


To solve log(x) = e/x, we first write the equation in the form f(x) = 0. This can be done by subtracting 
e/x from both sides of the equation, yielding log(x) — e/x = 0. We can rewrite this in the equivalent 
but slightly more convenient form x log(x) — e = 0. Since the left-hand side of this equation is un- 
defined for x <= 0, and we can’t guarantee that the search for a root will not venture into this region, 
we will consider instead the equation |x| log|x| — e = 0. This equation has exactly the same positive 
solution(s) as the first equation, but this equation makes sense for both positive and negative (but non- 
zero) numbers. The program below includes a user-defined function that computes the left-hand side of 
this equation, and uses FHROOT to find a solution of the equation. 





10 STD 


20 DEF FNF(X) This user-defined function computes the left- 
hand side of the equation. 


30 FNF =ABS(X)*LOG(ABS(X))—EXP(1) 

40 END DEF 

50 INPUT A,B These will be the initial guesses. 
60 R=FNROOT(A,B,FNF(FVAR)) 

70 DISP ‘R =';R 

80 DISP ‘FNF(R) =';FVALUE 

90 DISP ‘FGUESS=';FGUESS 


To use the program we must decide on initial guesses. Although the initial guesses need not be in 
increasing order, or even distinct, a choice of initial guesses that surround a root will produce results 
more quickly in general. Noting that if |F MAR |<1 then FHF «FAR» will be negative and if FAR is 
large (say, 100) then FHF ¢F\AR» will be positive, we can choose .5 and 100 for our initial guesses. 





Key in the program and [RUN] it, and when prompted with ? respond with . 5, 146 [END LINE), which 
supplies the initial guesses. The computer will then display 

















Since FHF «F® = &, the value given is an exact root for FHF. 
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Additional Information 


Choosing Initial Estimates 


When you use FHF OT to find roots of equations, the initial estimates determine where the search for 
a root will begin. If the two estimates surround an odd number of roots (signified by their function 
values having opposite signs), then FHF.OOT will find a root between the estimates quite rapidly. If the 
function values at the two estimates do not differ in sign, then FMF OMT must search for a region 
where a root lies. Selecting initial estimates as near a root as possible will tend to speed up this search. 
If you merely want to explore the behavior of the function near the initial etimates (such as to deter- 
mine if there are any roots or extreme points nearby), then specify any estimates you like. 


Another thing to consider is the range in which the equation is meaningful. In solving f(x) = 0, the 
variable x may only have a limited range in which it is conceptually meaningful as a solution. In this 
case, it is reasonable to choose initial estimates within this range. Frequently an equation that is ap- 
plicable to a real problem has, in addition to the desired solution, other roots that are physically 
meaningless. These usually occur because the equation being analyzed is appropriate only between 
certain limits of the variable. You should recognize this restriction and interpret the results 
accordingly. 


Interpreting Results 


FHROOT always evaluates the function at the value returned, as described above. This enables you to 
interpret the results. There are two possibilities: the value of the function at the value returned by 
FHROOT in close to 0; or the value of the function at the value returned by FHF OOT is not close to 0. 
It is up to you to decide how close is close enough to consider the value a root. 


If the function value is too large, then the information returned by the keyword FGLIES 5, together 
with information already considered, is sufficient to determine the general behavior of the function in 
the region. For example, suppose that FHF:C0T is used to find a root of a function—say, f(x)—and the 
value returned is r. If |F'\ALUE| is too large to consider r a root, then there are several possibilities. 








If F'\VALUE and f(F GUESS) have the same sign, then r is either an approximation to a local minimum 
of |f(x)| or in a region where the graph of f(x) is horizontal. 
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Case a Case b 


In the two cases above, FHF:O0T sees no tendency of f(x) to decrease in absolute value, and so to cross 
the x-axis. It will then try to approximate a local extreme point, if any. This approximation can be 
resolved to further precision by further executions of FMROOT, using r and FGUESS as initial es- 
timates. Repeated execution of FMROOT in this manner will converge to the extreme point in many 
cases. The idea is that FHROOT can be used to find local extreme points, or the information about 
where the extreme points are can be used to re-direct the search elsewhere, in hope of finding a root. 





When |F') FLUE | is too large to consider r a root, another possibility is that F/ALUE and f(F GUE 
have different signs. In this case it would appear that there is a root between, because for the function 
to change signs it should cross the x-axis. Typically, when F NF0I1T finds two guesses on opposite sides 
of the x-axis, it only stops after it has resolved them to two consecutive machine numbers. In this case 
there is no machine representable number between r and F GES. Thus, the behavior of the function 
cannot be determined between r and FGLIES 5. To interpret such results, you should be aware of these 


situations. 
er ao 
e pie 















| 
| 
| 
| 
| 
| 
| 
4 


Case 1 Case 2 Case 3 
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In case 1, r and FGU are the best approximations to the root that are representable on the ma- 
chine. Case 2 looks exactly the same to FHFOOT, but there is no root—there is a jump discontinuity 
instead. In case 3 there is a pole, which can look like a root if a guess on each side of the pole is found. 
FHROOT returns information in FGUESS and the F\\ALUE to help you isolate situations where 
convergence is to a pole. 








Decreasing Execution Time 


The exponent range of your HP-71 is +499 (except when TRAF «HF >» = 2, which effectively extends 
the negative exponent range to —510). This allows for sensitive observation of the behavior of a func- 
tion, even very close to a root. FHROOT takes advantage of this dynamic range by not accepting a 
guess as a root until the function value underflows, is zero, or two consecutive machine representable 
numbers that bracket a root are found. The cost of this precision is that, occasionally, it may take quite 
a while to obtain such precision. If this high degree of sensitivity is not required, then you may wish to 
set a smaller tolerance. For example, you may only need to know a place where the function is less than 
1E—20. This is accomplished in your function definition by checking the value of the function before 
assigning it to the function variable and setting the function variable to zero if the computed value is 
smaller than the desired tolerance. For example, suppose you wanted to find any root of f(x) = x‘, and 
|f(x)| <= 1E—32 is acceptable as a root. Here is a program you can use. 


10 STD 

20 DEF FNF(X) Multiline function definition of f(x) = x*. 

30 F=X%4 

40 IF F<=1.E—32 THEN FNF=0 ELSE Checks error tolerance and sets the function 

FNF=F value accordingly. 

50 END DEF 

60 DISP FNROOT(2,3,FNF(FVAR)) Computes and displays the root. 

70 DISP FVALUE Displays the function value at the root. 
Input/Result 
RUN 





In this example, if the tolerance technique were not used, execution would last much longer. This is 
because the computed function will not underflow until x is very small, since the root is at zero and the 


distribution of machine-representable numbers is very dense close to zero. So FMR OCT has a lot of 
guesses to try before finding one it can accept as a root. 
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An alternate approach to decreasing execution time is to translate the function so that the root is not 
so near zero, compute the root of the translated function, then translate the root back. This will de- 
crease the time to find roots of certain functions with roots close to zero, but will generally decrease 
the accuracy of the roots found. Here is a sample program for f(x) = x*. 


10 STD 

20 DEF FNF(X)=(X—1)*4 This is x* translated by 1. 

30 R=FNROOT(3,4,FNF(FVAR)) Computes the root. 

40 DISP R—-1 Translates the root back and displays the root 


and function value. 
50 DISP FVALUE 


Finally, there is a technique that may improve the speed and accuracy of FHF OOT. Any equation is 
typically one of an infinite family of equivalent equations with the same roots. However, some may be 
easier to solve than others. For example, the two equations f(x) = 0 and exp(f(x)) — 1 = 0 have the 
same real roots, but one is almost always easier to solve. When f(x) = x4 — 6x — 1, the first is easier; 
but when f(x) = In(x4 — 6x —1), then the second is easier. While FHROCT has been designed to 
provide accurate results for a wide range of problems, it is worthwhile to be aware of such possibilities. 


Suspending FHEGOT With 


If none of the arguments of FW 0T contain multi-line user-defined function calls, pressing will 
not terminate the action of FMFOOT until intermediate results are saved. In particular, FHRM0T will 
return and save the current F \/AR as though it were a root, it will save the previous guess as FGWESS, 
and it will save the value of f(x) at the current F\AF as the value of F\!ALUE. Only then will the 
action of F HR OOT stop. 








if, on the other hand, there are one or more multi-line user-defined functions as arguments for 
FNROOT (that is, if the calculation of FHROOT involves several BASIC program lines), [ATTN] will be 
ignored until a multi-line user-defined function is called. Execution will then halt at a Ine of one of the 
user-defined functions. This gives you the ability to examine relevent values, such as the current value 
of FAR, then continue the execution of FHF OCT if you wish. 





In addition, if there are multi-line user-defined functions as arguments for FMR COT, then fatal errors 
within the user-defined function do not destroy the FHROOT environment, giving you the exact same 
correct and continue capability as with any other HP-71 call to multi-line user-defined functions. 


CALC Mode 


You cannot execute FHEOOT directly or indirectly in CALC mode. For instance, suppose your current. 
file contains a single-line user-defined function F $F whose definition contains an FHF OO T keyword. If 
you attempt to execute FHF in CALC mode, an error will result. 
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Nesting Rules 


If the third argument F of FHFOOT defines a formula whose evaluation encounters another FMR O0T 
keyword, we say that the two FHFOOT keywords are nested. Up to five FHROOT keywords can be 
nested in this way. 


As an example of FHF: T nesting, consider the following program that solves f(x,y) = x2 + y? — 2x 
— 2y + 2 for x and y. 





10 STD 
20 DEF FNF(X,Y)=X*2+Y*2—2*X—2xY+2 Defines the function whose solution is sought. 
30 DEF FNG(X) Lines 30 through 60 define a one variable 
40 R=FNROOT(—4,4,FNF(X,FVAR)) function f(x) that receives a fixed x value (F VAR) 
50 FNG=FVALUE from line 70. 
60 END DEF 
70 DISP FNROOT(—3,3,FNG(FVAR));R If this FHROOT function receives a nonzero 
result from line 50, it selects another x value for 
the FHROOT in line 40 to try. If it receives a 
zero result, a solution for f(x,y) is found. 
Input/Result 
RUN 











The x and y values returned by the FMF OO T 
function in line 70. The x value is displayed on 
the left. 









The closest FFT comes to the true y value, one, is . 32 3, since these x and y values 
satisfy the objective of FHROOT. This objective is to find x and y values for which the computed value 
of f(x,y) is zero. 


F\ALUE (END LINE 


7 


These values for x and y when used in f(x,y) give 
& as the result. 
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A common use for FHF UOT nesting is determining minima. To demonstrate this application, we’ll 
modify the above function f(x,y) by adding one to the expression, thereby ensuring that the function 
has no solution, since the paraboloid represented by the modified function no longer intersects the xy 
plane. The only program modification is in line 20: 


20 DEF FNF(X,Y)=X*2+Y*2—2X—2*Y +3 


All other program lines are unchanged. 


The earlier nested FHROOT program required about 20 seconds to reach a solution. Since FHROOT 
takes special care to make sure a true minimum is found, the modified program requires about 312 
minutes to find and display the x and y values whose use in /(x,y) result in a function minimum. 


Input/Result 
RUN 


The x and y values that give a minimum for the 
Bidd44 modified function. 





FVALUE Displays the value given by the modified func- 


tion using these x and y values. 


There is no need to wait the full 34% minutes for a result. As explained on page 97, you can suspend an 
executing FMF COT function and then display interim results. If two consecutive inspections of interim 
results show insignificant change, you might wish to accept them as having satisfactory accuracy. 


Use of User-Defined Functions 


If the third argument of an FHF OMT function evaluates any user-defined function, then you must 
execute the FHFOOT function as a program statement, not from the keyboard. Also, if FHRCIOT is 
suspended while executing, you cannot execute a user-defined function from the keyboard, in either 
BASIC or CALC mode. 
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Numerical Integration 


Keywords 


You can use the keywords in this section to evaluate the integral of a function of from one to five 
variables between definite limits to an accuracy of your choosing. 


Throughout most of this section, the operation of these keywords will be described for a one-variable 
function. Multi-variable functions are covered under the topic Nesting Rules - Volume Integration, 


pages 109-110. 


The keyword IH TEGRAL can be used from the keyboard or inside a program to calculate the integral 
of the function, provided the keyboard line or program contains the function definition. 


The keywords IBGUHO and 1\/ALUE give you additional flexibility in the evaluation of the integrals. 
INTEGRAL, ITBOUND, and I\!ALUE are numeric-valued, so they can be used alone or in combination 
with other functions and variables to form numeric expressions. A fourth keyword, I! AF, represents 
the variable (or one of the variables) of integration in the function being integrated by INTEGRAL. It 
also contains the most recent sampling point used by an executing INTEGRAL. 


INTEGRAL Definite Integral 
IHTEGRALSA,B,E,F? 
where A, B, E, and F are real numeric expressions. 


Returns an approximation to the integral from A to B of F. The relative error E (rounded to the range 1E- 
12<=E<=1) indicates the accuracy of F and is used to calculate the acceptable error in the 


approximation of the integral. 


This integral approximation can be: 
e An approximation to the integral that is as accurate as the relative error E would allow. 
© The last of 16 approximations to the integral, which have sampled the integrand at 65535 points 
without meeting the convergence criterion. 
© The best current approximation to the integral returned when [ATTN] is pressed and when F does 
not call a multi-line user-defined function. 
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INTEGRAL (continued) 

IHTEGRAL generates a sequence of increasingly accurate approximations to the definite integral. If 

three successive approximations are within the acceptable error of each other—the first is close to the 

second and the second is close to the third—they are considered to have converged and the third 

approximation is returned as the value of the definite integral. If a total of 16 approximations are cal- 
culated without converging, the sixteenth is returned. 


Not usable in CALC mode. Refer to page 111 for more information about IMTEGRAL and CALC mode. 


Refer to pages 109-111 for information about IHTEGRAL nesting (volume integration) and about the: 
interactions between IHMTEGRAL and and between IMTEGRAL and user-defined functions. 


IVAR Integration Variable 


IVAR 


Represents the variable of integration in the formula defining F, the last argument of [MTEGRAL. 


| Also contains the most recent sampling point used by a running IMTEGRAL. 





Can be used in CALC mode. 


IVALUE Last Result of INTEGRAL 
IWALUE 


Returns the last approximation computed by the 1MTEGRAL keyword. If the [ATTN] key was pressed 
or the operation of IHTEGRAL was otherwise interrupted, then I\/ALUE returns the value of the 
current approximation to the integral. Otherwise, 1\/HLUE returns the same value that INTEGRAL 
| last returned. 


IVALUE retains its value (even if your HP-71 is turned off) until another IHTEGREAL is computed. 











| 
Can be used in CALC mode. 
l 
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IBOUND Error Approximation for INTEGRAL 





TBOUHD 





Returns the final abso/ute error estimate for the definite integral most recently computed by 
INTEGRAL. 


e A positive value for 1 B00 means that the approximations converged. 


e A negative value for TE 0 LIND means that the approximations didn’t converge, so that the value 
returned by INTEGRAL may not be representative of the true value. 





Like ['\JALUE, IBC UND retains its value (even if the HP-71 is turned off) until another THTEGRAL. is 
computed. Unlike I'ALUE, the value of [OUND has no relation to the current approximation to the 
integral if the operation of INTEGRAL is interrupted. 


Can be used in CALC mode. 





To integrate a function between bounds, you can follow these steps: 


1. Write down an expression that represents the function to be integrated. 


. Substitute the keyword I \A everywhere in the expression for the variable of integration. 


ew 8 


. Use this expression as the fourth argument F of INTEGRAL. 


4, Use the lower and upper bounds of integration as the first and second arguments A and B of 
INTEGRAL, respectively. 


o 


. Choose a value for the third argument F of INTEGRAL that represents an estimate of the relative 
error in the computation of the integrand. Any value for E is always rounded to the range 
[1E—12,1]. Thus, E should satisfy, after rounding 

|TRUE INTEGRAND — COMPUTED INTEGRAND| 3S 
|COMPUTED INTEGRAND| 








Since INTEGRAL has no way of knowing what the true value of the function is intended to be, 
only you can supply this estimate. For many purely mathematical functions (STH, EXF, poly- 
nomials, etc.) and modest limits of integration, full 12 digit accuracy can be returned so that a 
value for E around 1E-12 should be suitable. 
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The operation of INTEGRAL and IBQUNO can be described more precisely as follows. 


1. Based on a relative error of E for the specified function, the computer calculates an error tolerance 
for the integral of the specified function. If f(X) is the “true” function that F approximates, then 
choose E£ such that 


IF — f(X)| . 
—— =F 
| 


for all X in the interval of integration. Your input for E is rounded to the range 1E—12<E<1, 


For example, if F is derived from experimental data with N significant digits, let E equal 10—%. 


2. The computer calculates a sequence of approximations /;, to the integral of the specified function. 
The difference between successive approximations is compared to the error tolerance for the 
integral. 


3. A value for the integral is returned when 


e The approximations J, have converged. Convergence is determined using J, defined as the kth 
approximation to the integral of E+|F| over the same interval of integration. J), represents the 
error inherent in the computation of /;,. 


The approximations J, are judged to have converged to I,, if 
[e — Th-al<Ie 


fork = n — landk = n. The value of J, is then returned by INTEGRAL; a positive value for 
the error estimate will be returned by IBOQUND. 
or when 


e The computer has evaluated J, through J; but the convergence criterion is still not met. I), is 
then returned by INTEGRAL; a negative value for the error estimate will be returned by 
IBOUND. 








Examples 
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Integrating f(x) = x2—2 (IHTEGRAL, IVA Ft) 


The following six examples illustrate various ways I HTEGRAL and IA can be used to integrate the 
function x? — 2 from i to 2. The first and sixth examples show the solution. 


Example One. 


Input/Result 
THTEGRAL“®1,2,1E-11,I1VAR 


END LINE 











Example Two. 


10 DISP INTEGRAL(COS(0),LOG2(4), 
1E—11,IVAR*2—2) 


Example Three. 


10 DEF FNG=IVAR*2—2 
20 DISP INTEGRAL(1,2,1E—11,FNG) 


Example Four. 


10 DEF FNF(X)=x*2—2 
20 DISP INTEGRAL(1,2,1E—11,FNF(IVAR)) 


Example Five. 


10 DEF FNH 

20 FNH=IVAR*2—2 

30 END DEF 

40 DISP INTEGRAL(1,2,1E—11,FNH) 


IHTEGRAL can be used from the keyboard as 
well as in a program. 


The limits of integration can be expressions. 


The fourth argument of IMTEGRAL can be an 
expression or a reference to a user-defined 
function. 


I\!AF can appear in the user-defined function, 
as above, or in the fourth argument of 
IHTEGRAL. 


The user-defined function can consist of one or 
several lines. 





106 Section 11: Numerical Integration 


Example Six. 


10 DEF FNJ(X) 
20 FNJ=X*2—2 


30 END DEF 

40 DEF FNF (X)=2*x 

50 DISP INTEGRAL(1,FNF(1),1E—11, The first, second or third arguments of 
FNJ(IVAR)) INTEGRAL can be references to user-defined 


functions. 
60 DISP IBOUND 


Input/Result 
RUN 





The resulting integral. 


The absolute error estimate for the resulting 
integral. Since it’s positive, the approximations 
converged. 





Integrating f(x) = e*—2 (IHTEGRAL, IVAR, IVALLE) 


This example features I''ALUE. This function returns the most recent integration approximation and 
is updated even while the execution of IMTEGRAL is in progress. After the execution of IMTEGRAL 
is completed, I\ALUE returns the same value returned by INTEGRAL. 


You can watch the progress of integral approximations by displaying | \!ALUE during the execution of 
INTEGRAL. This is demonstrated by the following program, which integrates the function e“—2 from 
one to three. The error bound used is 1E-12. 


10 Y=IVALUE ‘’ = value of I\/ALLIE when program starts 
(assumes I'\JALUES is set from a previous 
INTEGRAL). 


20 DEF FNF(X) 
30 IF IVALUE=Y THEN 50 Displays I!/ALUE only if it has changed. 
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40 DISP IVALUE @ Y=IVALUE 
50 FNF=EXP(x)—2 
60 END DEF 


70 DISP INTEGRAL(1,3,.000000000001, 
FNF(IVAR)) 


Input/Result 


RUN 





First displayed value of I\/ALUE. 


Last displayed value of I\/ALLIE. 
Value of INTEGRAL. 





Integrating f(x) = exp(x9+4x2+x+1) (INTEGRAL, IVAR, TBOUND, 
IVALUE) 


To find the integral from 0 to 1 of the function 
f(x) = exp(x3 + 4x2 + x + 1) 


you can use the following program. 


10 DEF FNF(X)=EXP(X*3 + 4*X*2+ X+1) The user-defined function FHF. 


20 INPUT E Gets the relative error we expect in FHF as 
compared with /. 


30 DISP ‘Integrating’ 

40 X=INTEGRAL(0,1,E,FNF(IVAR)) 
50 BEEP 

60 DISP ‘Integral =’;X 

70 DISP ‘The approx. error =" 

80 DISP IBOUND 
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After you key in the program, run it using the following keystrokes. 

















Input/Result 
RUN 
? & The prompt to enter the relative error of the 
function. 
LE~5 (END LINE Although our function is accurate to one part in 











10!2, we can say that it is less accurate (in this 
case, one part in 10°) so that the computation 
will finish more quickly. 


The value of the integral is 104.2911 + (1.04 x 
10-8). 








I WALUE (END LINE 


IVALUE gives the value of the last computed 
integral. 





Integrating C(T) = a + bT (INTEGRAL, IVAR, IBQUND) 


You can use 1TEGRAL to compute the amount of heat required to heat one gram of gas at a constant 
volume from one temperature to another. The amount of heat needed, Q, is given by the formula 


T: 
a= fran ar 


where C(T) is the specific heat of the gas as a function of temperature, 7'l is the starting temperature, 
and 72 is the final temperature. 


If C(T) = a + bT, where a and b are experimentally determined to be a = 1.023E~2 and b = 2.384E 2 
with four significant digits, then we can compute the relative error of C(T) to be approximately 5E—4. 
The program below prompts you for the initial and final temperature in degrees Kelvin and then com- 
putes the heat needed to raise the temperature of the gas from the initial to the final temperature. 
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10 DEF FNC(T)=.01023+ .02384*T The user-defined function that calculates the 
specific heat. 


20 INPUT ‘Initial, final T (K)?’;T1,T2 

30 DISP ‘Integrating’ 

40 Q=INTEGRAL(T1,T2,.0005,FNC(IVAR)) Computes the integral. 

50 DISP ‘Heat needed =’;Q;‘+ —’;IBOUND Displays the answer and the approximate error. 


To find the heat needed to raise the temperature from 300°K to 310°K, type in the program and use the 
following keystrokes. 


Input/Result 
RUN 


Additional Information 


Nesting Rules—Volume Integration 


If the fourth argument F of INTEGRAL defines a formula whose evaluation encounters another 
INTEGRAL keyword, we say the two INTEGRAL keywords are nested. Up to five INTEGRAL 
keywords can be nested in this way. A program that nests two INTEGRAL keywords can determine 
volumes. 


As an example of IMTEGRAL nesting, consider the following program that integrates f(x,y), where 
f(x,y) = x? + 2y, over the square 0 <x <1, 0 <y <1. That is, the program evaluates 


I : e “x, y)dyde. 
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10 DEF FNF(X,Y)=X*2+2*Y Defines the function whose integral is sought. 

20 DEF FNG(X)=INTEGRAL(0,1,1E—6, For each value of ™, integrates a slice parallel to 
FNF(X,IVAR)) the y-axis. 

30 INTEGRAL(0,1,1E—6,FNG(IVAR)) Sums all of the contributions from the slices 


parallel to the y-axis. 


Input/Result 
RUN 














The volume returned by the INTEGRAL fune- 
tion in line 30. 





The answer is exact even though TE OUNMO only 
predicts six correct digits. 





The following example demonstrates the use of INTEGRAL to evaluate the integral 
r/2ry . 
J J sin(x) dxdy 


Input/Result 
RAOQIAWS LEND LINE 


IHTEGRAL¢@, PI TE-3, Note that the first 1\/AF is the integration vari- 
INTEGRA » IVAR; 1E=3, able of the outside INTEGRAL. 


> (END LINE) 


























The true answer is /2 — 1 (approximately 
-5707963268). 
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Suspending IHTEGRAL With 


If none of the arguments of INTEGRAL contain multi-line user-defined function calls, pressing [ATTN 
will not terminate the action of INTEGRAL until intermediate values are saved. In particular, 
INTEGRAL will save and return the current I'ALUE as though it were the integral, and will make 
negative the current value of IEGUHO. Only then will the action of INTEGRAL stop. 














If, on the other hand, there are one or more multi-line user-defined functions as arguments for 
INTEGRAL (that is, if the calculation of IHTEGRAL involves several BASIC program lines), [ATTN 
will be ignored until a multi-line user-defined function is called. Execution will then halt at a line of 
one of the user-defined functions. This gives you the ability to examine relevent values, such as the 
current value of I \!ALUE, then continue the execution of INTEGRAL if you wish. 














In addition, if there are multi-line user-defined functions as arguments for INTEGRAL, then fatal 
errors within the user-defined function do not destroy the INTEGRAL environment, giving you the 
exact same correct and continue capability as with any other HP-71 call to multi-line user-defined 
functions. 


CALC Mode 


You cannot execute IMTEGRAL directly or indirectly in CALC mode. For instance, suppose your cur- 
rent file contains a single-line user-defined function FHF whose definition contains an INTEGRAL 
keyword. If you attempt to execute FHF from CALC mode, an error will result. 


Use of User-Defined Functions 


If the fourth argument of an INTEGRAL function evaluates any user-defined function, then you must 
execute the INTEGRAL function as a program statement, not from the keyboard. Also, if IMTEGRAL 
is suspended while executing, you cannot execute a user-defined function from the keyboard, in either 
BASIC or CALC mode. 


Overview of Numerical Integration 


Numerical integration schemes generally sample the function to be integrated at a number of points in 
the interval of integration. The calculated integral is simply a weighted average of the function values 
at these sample points. Since a definite integral is really an average value of a function over an infinite 
number of points, numerical integration can produce accurate results only when the points sampled are 
truly representative of the function’s behavior. 








If the sample points are close together and the function does not change rapidly between two consec- 
utive sample points, then the numerical integration will give reliable results. On the other hand, 
numerical integration will not produce good answers on a function whose values vary wildly over a 
domain that is small in comparison with the region of integration. Other errors that can affect the 
result of a numerical integration include the round-off errors typical of any floating point computation 
and errors in the procedure that computes the function to be integrated. 


Handling Numerical Error 


The INTEGRAL keyword requires specification of an error tolerance E for its operation. This error 
tolerance is taken to be the relative error of the computed function as compared with the “true” func- 
tion to be integrated. The error tolerance is used to define a ribbon around the computed function and 
the “true” function should then lie inside this ribbon. If the “true” function is f(x) and the computed 
function is F(x), then 


F(x) — Error (x) < f(x) < F(x) + Error (x) 


where Error (x) is half the width of the ribbon at x. 


f(x) 





We can then conclude that 


£ f(x) dx = fe F(x) dx + § Error (x) dx 


where the third integral is just half the area of the ribbon—that is, integrating the computed function 
instead of the actual function can introduce an error no greater than half the area of the ribbon. 
INTEGRAL estimates this error while computing the integral; IE OUND gives you access to the 
estimate. 





eS ms 
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Choosing the Error Tolerance 


The accuracy of the computed function depends on three factors: 


. The accuracy of empirical constants in the function. 


, The degree to which the function may accurately describe a physical situation. 
, The round-off error introduced when the function is computed. 


Functions like cos(x — sinx) are purely mathematical functions. This means that the functions contain 
no empirical constants, and neither the variables nor the limits of integration represent any actual 
physical quantities. For such functions you can specify as small an error tolerance as desired, provided 
the function is calculated within that error tolerance (despite round-off) by the BASIC function. Of 
course, due to the trade-off between accuracy and computation time, you may choose not to specify the 
smallest possible error tolerance. Any specified error tolerance is rounded to the range [1E—12, 1]. 


When the integrand relates to an actual physical situation, there are additional considerations. In these 
cases, you must ask yourself whether the accuracy you would like in the computed integral is justified 
by the accuracy of the integrand. For example, if the function contains empirical constants that 
approximate the actual constants to three digits, then it may not make sense to specify an error tol- 
erance smaller than 1E—3. 


An equally important consideration, however, is that nearly every function relating to a physical situa- 
tion is inherently inaccurate because it is only a mathematical model of an actual process or event. A 
mathematical model is typically an approximation that ignores the effects of factors judged to be in- 
significant in comparison with the factors in the model. 


For example, the equation s = s’ — (.5)gt?, which gives the height s of a falling body when dropped 
from an initial height s’, ignores the variation with altitude of g, the acceleration due to gravity. 
Mathematical descriptions of the physical world can provide results of only limited accuracy. If you 
calculated an integral with an accuracy greater than your model can support, then you would not be 
justified in using the calculated value to its full (apparent) accuracy. It makes sense to supply an error 
tolerance that reflects any inaccuracies in the function, or the INTEGRAL keyword will waste time 
computing to a level of accuracy that may be meaningless. Further, the value returned by IBOUNMO 
may not be significant. 


If f(x) is a function relating to a physical situation, its inaccuracy due to round-off is typically very 
small compared to the inaccuracy in modelling the situation. If f(x) is a purely mathematical function, 
then its accuracy is limited only by round-off error. Precisely determining the relative error in the 
computation of such a function generally requires a complicated analysis. In practice, its effects are 
determined through experience rather than analysis. 








Handling Difficult Integrals 


Integrating on Subintervals. A function whose values change substantially with small changes in 
its argument will likely require many more points than one whose values change only slightly. This is 
because the behavior of the function must be adequately represented by the sampling. If a function is 
changing more rapidly in some subintervals of the interval of integration than in others, you can sub- 
divide the interval and integrate the function separately on the smaller intervals. Then the integral 
over the whole interval is the sum of the integrals over all the subintervals, and the error of the integral 
is the sum of the errors of the integrals over the subintervals. 


The algorithm used by INTEGRAL makes a reasonable decision during execution of how many points 
to sample, based on the behavior of the specified integrand on a particular interval. When the interval 
of integration is split up, each subinterval can be handled according to the function behavior on that 
subinterval alone. This results in greater speed and precision. 


For example, to integrate f(x) = (x? + 1E—12)” from x = —3 to x = 5 using an error tolerance of 
1E—12, it speeds up execution to subdivide the interval at x = 0, where f(x) has a sharp bend in its 
graph. Because f(x) is very smooth on the subintervals (—3, 0) and (0, 5), the integrals over these 
subintervals can be evaluated quickly. 


5 0 5 
Jo, fe) dx = fo, fx) de + [fo dx 
The following program computes this integral on the two subintervals and then combines the results. 


10 DEF FNF(X)=SQR(X*X+ 1E—12) We will use *#* rather than } because Hq 
is more accurate. An analogous situation generally 
occurs for any integer power of a variable. 


20 |I=INTEGRAL(—3,0,1E—12,FNF(IVAR)) Integrate over the first subinterval. 

30 E=IBOUND Save the error to add in later. 

40 DISP “Integral ="; 

50 DISP |+INTEGRAL(0,5,1E—12,FNF(IVAR)) The sum of the first and second integrals. 


60 DISP “Error =";E + IBOUND Compute the relative error by adding the two 
errors together. 





You can run this program by keying it in and then pressing [RUN]. The following will then appear in 
the display. 
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When the interval is subdivided, I HT EGRAL computes the answer in a few seconds. Without subdivid- 
ing the interval, execution may take a long time. 


Subdividing the interval of integration is also useful for functions with a singularity in the interval. 
The singularity may consist of one or more points where the function is undefined or has a sharp 
corner point. 


For example, the integral 


2 dx be, L dx 2 dx 
ic Gi may be split into } GD + f @— Db 


to avoid evaluating the function at x = 1, where it is undefined. You can now integrate the function on 
each subinterval because x = 1 is an endpoint of each subinterval, and IMTEGRAL does not sample at 
an endpoint. 


Similarly, the function \/|x — 1|, has a sharp corner point at x = 1. 





Suppose you need to integrate this function from 0 to 2. You can increase the speed and accuracy of the 
computation by integrating separately on the subintervals (0, 1) and (1, 2), because the function is 
smooth on each of these subintervals. 


Transformation of Variables. A second method of handling difficult integrands is by transforming 
the variable. When the variable in a definite integral is transformed, the resulting definite integral may 
be easier to compute numerically. Consider the integral 


i = = ty) 


The derivative of the integrand approaches infinity as x approaches 0, as shown on the left below. The 
substitution x = u? stretches the x-axis and causes the function to be better behaved, as shown on the 
right. 
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0.1 


2u? u 
(u + 1)\(u —1) Inu 


0.1 





You can now evaluate the integral that results from this substitution: 
1 2 
(ia 
0 \(u + 1)(u — 1) Inu | 


(Do not replace (u + 1)(u — 1) with u? — 1; as u approaches 1, u2 — 1 loses half of its digits to 
roundoff, yielding a final result that is too large.) 





As a second case requiring substitution, consider the following function. Its graph has a long tail 
stretching much farther than the main body (where most of the area is). 
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Although a very thin tail may be truncated without greatly degrading accuracy, this function has too 
wide a tail to ignore when calculating 
ete 
-t x +197) 


if is large. In general, the compressing substitution x = b tan u maps the entire real line into (— 71/2, 
7/2) and maps subsets of the real line into subsets of (—7/2, /2). For b = 1E—5 the substitution 
becomes x = 1E—5 tan u and the integral becomes 


5 tan '(t/b) 
a0 tee 


which is readily computed for very large t. 


This compressing substitution is also a standard way to deal with infinite intervals. For example, 


es de ag Re 
te 2 +1977 1° ate a 


Tn some cases the tail can be chopped off. Consider the function exp (—x2). This functions underflows 
(that is, gives a result of zero in machine arithmetic) for x > 34. Thus, 


& e* dx ® tig e “dx 


Therefore, when dealing with infinite integrals you can cut off the tail if it is insignificant, but you 
should use a compressing substitution if it is not. 
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About the Algorithm 


The Math Pac uses a Romberg method for accumulating the value of an integral. Several refinements 
make it more effective. Instead of equally spaced samples, which can introduce a kind of resonance or 
aliasing that produces misleading results when the integrand is periodic, IHTEGF AL uses samples 
that are spaced nonuniformly. Their spacing can be demonstrated by substituting 


es ie 
x= Suu into f° fede 


and then spacing u uniformly. Besides suppressing resonance, the substitution has two additional bene- 
fits. First, no sample need be taken from either endpoint of the interval of integration unless the inter- 
val is so small that points in the interval round to an endpoint. As a result, an integral like 


1 sin x 
if} ae 
0 x 


will not be interrupted by division by zero at an endpoint. Second, IHTEGF AL can integrate functions 
whose slope is infinite at an endpoint. Such functions are encountered when calculating the area en- 
closed by a smooth closed curve like x2 + f2(x) = R. 


In addition, INTEGRAL uses extended precision. Internally, sums are accumulated in 15-digit num- 
bers. This allows thousands of samples to be accumulated, if necessary, without losing any more signifi- 
cance to round-off than is lost within your function. 


During the computation, INTEGRAL generates a sequence of iterates that are increasingly accurate 
estimates of the actual value of the integral. It also estimates the width of the error ribbon at each 
iterate. IHTEGRAL stops only after three successive iterates are within the computed error of each 
other or after 16 iterations have been performed without this criterion being met. 


In the latter case the function will have been sampled at 65,535 points. The value returned by 
I BOUND will be the negative of the computed error to signify that the returned value of the 
INTEGRAL is likely not within the error tolerance of the actual value. Typically, you should then split 
up the interval of integration into smaller subintervals and integrate the function over each of the 
subintervals. The integral over the original interval will then be the sum of the integrals over the 
subintervals. In this way, up to 65,535 points can be sampled on each subinterval, thus computing the 
integral to greater precision. 





In summary, IHTEGRAL has been designed to return reliable results rapidly and in a convenient, 
easy-to-use fashion. The above theoretical considerations discuss problems with numerical integration 
in general. The [HTEGRAL keyword is capable of handling even difficult integrals with relative ease. 
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Finding Roots of Polynomials 
















eyword 


keyword in this section finds all solutions—both real and complex—of P(x) = 0, where P is a 

lynomial of your choice with real coefficients. If P is a polynomial of degree n there will be n (not 
necessarily distinct) solutions of this equation, so this keyword resembles an array operation in its 
rmat. 


‘To use this keyword to find the solutions of the equation P(x) = 0, where 
P(x) = apx"” + Qy_yx™-! + ... + ax + ag 


first store the coefficients a,, a, ;, ..., 4g in a real-type array with n + 1 elements in all. They should 
be stored in the order indicated above, with the coefficient of the highest power first and the constant 
term last. Aside from the total number of elements in the array, which indicates to the Math Pac the 
degree of the polynomial, the dimensions of the array are irrelevant. For example, the arrays 


65 4 
[6, 5, 4, 3, 2, 1], ,|4 3], and 
21 


re No Oo Fk TD OD 


all can represent the fifth degree polynomial 
6x5 + 5x4 + 4x3 + 3x2 + 2x + 1. 


The array in which you wish the resulting roots to be stored must be complex type to accept complex 
roots. If the polynomial whose roots are sought has degree N, and if the result array is a vector, it will 
be redimensioned to have N elements. If the result array is a matrix, it will be redimensioned to have 
N rows and one column. 


The degree of the polynomial whose roots you can find is limited only by the amount of memory you 
have available. 
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PROOT Roots of a Polynomial 





MAT R=PROOTCP> 


where P is a real vector or matrix with N + 1 elements, where N = degree of polynomial whose roots 
are sought, and R is a complex vector or matrix. 





If R is a vector, implicitly redimensions R to have N elements. If R is a matrix, implicitly redimensions R to 
have N rows and one column. R will be assigned the (complex) values of the solutions of the equation 
P(x) = 0 (where P is the polynomial of degree N whose coefficients are the values of the elements of P). 





To halt operation, press [ATTN] twice. 


Not usable in CALC mode. 





Example 
The following example finds all roots of the polynomial 


5Z® — 45Z° + 225Z4 — 425Z° + 170Z? + 370Z — 500 








OPTION BASE 1 [ENDLINE 
OIM ACPO Creates real vector for coefficients. 
COMPLEX Be 18> (ENDLINE Creates complex vector for roots. 


MAT INPUT A LEND LINE 





5,495,225 ,-425,;178, 370, =500 





END LINE 
MAT B=PROOTCA> (ENDLINE First redimensions the vector B to have six ele- 


ments, just large enough to contain the six (com- 
plex) roots of the six-degree polynomial. Then 
finds all roots and stores them in B. 
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HAT OS 











E [END LINE Displays all roots. 








Additional Information 


There are several methods of gauging the accuracy of the calculated roots. The first method is to cal- 
culate the value of the polynomial at the alleged root, and compare this value with zero. Although quite 
straightforward in theory, this has a number of drawbacks in practice. It may easily happen that the 
root calculated is the closest machine-representable number to a true root, but because the polynomial 
has such a large value for its derivative at this root, the value of the polynomial at the calculated root is 
very large. A simple example of this phenomenon is given by the polynomial 1E20x2 — 2E20. A true 
root is \/2; a calculated root is 1.41421356237, which is the machine-representable number closest to 
\/2. However, the value of the polynomial at this approximation to the square root of 2 is 
—1,000,000,000, a number that seems very far from zero. 


Another drawback of the above method is that because of the limited precision available in any numeri- 
cal calculation, the roundoff errors that occur in the calculation of the polynomial value may com- 
pletely eliminate the significance of the difference between the calculated value and zero. This is 
especially true when the polynomial is of large degree, has coefficients widely varying in size, or has 
roots of high multiplicity. 


A second method of gauging the accuracy of the calculated roots is to attempt to reconstruct the poly- 
nomial from these roots. If the reconstructed polynomial closely resembles the original, the roots are 
then judged to be accurate. This technique is less sensitive to the problems that affect the polynomial 
evaluation method. Of course, this method does not give information on the accuracy of an individual 
root. 


The program below asks you for a polynomial and calculates the roots of that polynomial using the 
PROUT keyword. If you wish, the program will also calculate the reconstructed coefficients from the 
c..culated roots. In addition, if desired, the program will compute the value of the polynomial at either 
a calculated root or any other real or complex value. 


Lines 10 through 200 drive the program and use the F'F (iT function to calculate the roots of the given 
polynomial. Lines 210 through 250 comprise the subroutine that evaluates the polynomial at any real or 
complex point. Horner’s method is used. 
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Lines 260 through 410 comprise the subroutine that reconstructs the coefficients from the calculated 
roots, Starting with the polynomial 1, the subroutine successively multiplies the polynomial by the 
linear factors (Z — R), where R is a calculated real root, or by the quadratic Z2 — 2REPT(R) + 
ABS(R)?, where R is a calculated complex root. (Note that CONJ(R) will also be a calculated root). 


10 OPTION BASE 0 @ INTEGER D,E 
@ DIM U$[4] @ DELAY 1 @ WIDTH 96 


20 INPUT “DEGREE? ";D 
30 DIM P(D),C(D) @ COMPLEX R(D—1) 


40 DISP “ENTER COEFFICIENTS " 
@ MAT INPUT P 


50 DISP “WORKING...” 
60 MAT R=PROOT(P) 


70 DISP “THE ROOTS ARE” @ DELAY 8 @ 
MAT DISP R @ DELAY 1 


80 U$=KEY$ @ INPUT 
“RECONSTRUCT? (Y/N) ”;U$ 


90 IF UPRC$(U$)="Y" THEN GOSUB 260 
ELSE 110 


100 DISP “RCON COEFFICIENTS ARE” @ 
DELAY 8 @ MAT DISP C @ DELAY 1 


110 US=KEY$ @ INPUT 
“EVALUATION? (Y/N) ";US 
120 IF UPRC$(US)#*Y" THEN 190 

ELSE COMPLEX Z 


130 INPUT “AT A ROOT? (Y/N) ";U$ 


140 IF UPRC$(U$)#*°Y" THEN INPUT 
"VALUE? ";Z @ GOTO 160 


OD is the degree of the polynomial. 


Array P will contain the coefficients of the 
polynomial in the order given previously, array R 
will contain the calculated roots, and array C 
will contain the reconstructed coefficients. 


Enter the coefficients. The leading coefficient 
should be nonzero for the program to work 
properly. 


Calculates the roots and stores them in array R. 


Displays the calculated roots. To continue the 
program after each root is displayed, press 


[END LINE}. 


If you wish, the program will reconstruct the 
coefficients from the calculated roots. 


The subroutine starting at line 260 performs the 
reconstruction and stores the reconstructed 
coefficients in array C. 

Displays the reconstructed coefficients. To 
continue the program after each display, press 
[END LINE}. 

If you wish, the program will evaluate the 
polynomial at a root or at any other point. 


The complex variable 2 will hold the polynomial 
value. 


The point may be either real or complex. 











150 DISP USING '#,"WHICH ROOT 
(1...",K,")”;D @ INPUT E 
@ Z=R(E-1) 

160 GOSUB 210 @ DISP “POLYNOMIAL 
VALUE IS” @ DELAY 8 @ DISP Zz @ 
DELAY 1 


170 US$=KEY$ @ INPUT 
“ANOTHER VALUE? (Y/N) ";U$ 


180 IF UPRC$(U$)=“Y" THEN 130 
190 INPUT “ANOTHER POLY? (Y/N) ";U$ 


200 IF UPRC$(U$)="Y" THEN 20 ELSE STOP 
210 COMPLEX B @ B=P(0) 


220 FOR K=1 TOD 

230 B=P(K)+Z*B 

240 NEXT K 

250 Z=B @ DESTROY B @ RETURN 
260 DISP “WORKING...” 


270 MAT C=ZER @ C(D)=1 

280 FOR L=1 TOD 

290 IF IMPT(R(L—1))#0 THEN 340 

300 FOR K=D—L TO D—1 

310 C(K)=C(K +1)—C(K)*REPT(R(L—1)) 

320 NEXT K 

330 C(D)= —C(D)KREPT(R(L—1)) @ GOTO 
400 


340 REAL B @ B=REPT(A(L—1))*2 
+-IMPT(R(L—1))*2 
350 FOR K=D-L~—1 TO D-2 
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Input the number of the root where you want 
the polynomial evaluated. 


The subroutine beginning at line 210 evaluates 
the polynomial at the given point or root. This 
value is then displayed. To continue, press 


[END LINE]. 


The program will evaluate the polynomial again 
if you wish. 


You can choose to start the program over again 
with a new polynomial. 


The polynomial evaluation subroutine uses 
Horner’s method. 


This line begins the coefficient reconstruction 
subroutine. Some rounding error may accumulate: 
during reconstruction, so even if the roots are 
exact, the reconstructed coefficients may not 
exactly coincide with the original coefficients. 
Creates polynomial 1 in array C. 


We use each calculated root in turn. 


Lines 300 through 330 multiply the current 
reconstructed polynomial by a linear factor. 


Lines 340 through 390 multiply the current 
reconstructed polynomial by a quadratic factor. 
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360 C(K)=C(K + 2)—2*REPT(R(L—1)) 
*C(K+1)+B*C(K) 

370 NEXT K 

380 C(D—1)=—2*REPT(R(L—1))*C(K +1) 
+B*C(K) 


390 C(D)=B*C(D) @ L=L+1 


400 NEXT L 


410 MAT C=(P(0))kC @ DESTROY B 
@ RETURN 


We increment L since we multiplied the 
polynomial by both the complex root and its 
complex conjugate. 


The reconstructed polynomial has leading 
coefficient 1 and so must be adjusted if the 
original leading coefficient was not 1. 


Example. If we wanted to find and evaluate the roots of the polynomial 


O+ x + xt + 8+ 2 + x41, 


we would run the program using the following keystrokes. 


Input/Result 
RUN 







DEGREE? @ 


© LEND LINE 













EHTER COEFFICIENTS 








1,9, 45 1ed5 4-1 [ENDLINE 
















HORE IMG... 
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The display scrolls to display the imaginary part 
of the first root. 





The real part of the first root. 


The imaginary part of the second root. 




















The real part of the second root. 








Fourth root: © 


Fifth root: « . 





Sixth root: ©. 
After the last root is displayed, continue the program by pressing [END LINE]. 


Input/Result 


RECOMSTRUCT? ¢yeHo £ Any response but * or ‘: is interpreted as “no.” 





‘Y [END LINE 
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ROOM COEFFICIENTS ARE 





isc atari : iid . a 


END LINE 


The coefficient of the x° term. 





Display the remaining five coefficients in the same way. These displayed coefficients are: 
Coefficient of x* term: i 


Coefficient of x° term: . 3339 





Coefficient of x? term: 1 
Coefficient of x! term: . 33399599995 
Coefficient of x° term: 1 


After the last coefficient is displayed, continue the program by pressing [END LINE]. 





Input/Result 


EVALUATION? ¢YeHo # 





‘Y [END LINE 








*' [END LINE 


WHICH ROOT 
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1 [END LINE 


POLYHOMIAL WALUE I5 





END LINE Continues the program. 


ANOTHER VALUE? <YeN> @ 





Y END LINE] 


AT A ROOT? {YeN> © 





MN [END LINE J 


€-,2,.95 (ENDLINE 


POLYMOMNIAL VALUE IS 





ANOTHER VALUE? ¢Y-N> 





M [END LINE 


AHOTHER POLY? ¢Y-H> @f 








4 [END LINE Ends the program. 
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About the Algorithm 


The Math Pac finds the roots of polynomials using Laguerre’s method, which is an iterative process. 
The Laguerre step at the iterate Z, for the polynomial P(Z) of degree N is 


—NP(Z,) 
P(Z,) + ((N — 1)? (P(Z,))? — N WN — 1) P(Z,) P"(Z,)1" 


The sign in the denominator is chosen to give the Laguerre step of smaller magnitude. Polynomials or 
their quotients of degree < 3 are solved using the quadratic formula or linear factorization. 


Laguerre’s method is cubically convergent to isolated zeros and linearly convergent to zeros of mul- 
tiplicity greater than one. 


The FOOT function is global in the sense that the user is not required to supply either an initial guess 
or a stopping criterion; in other words, no prior knowledge of the location of the roots is assumed. The 
FOOT function always attempts to begin its search (iteration) at the origin of the complex plane. An 
annulus in the plane known to contain the smallest magnitude root of the current (original or quotient) 
polynomial is constructed about the origin (using five theoretical bounds) and the initial Laguerre step 
is rejected if it exceeds the upper limit of this annulus. In this case, a spiral search from the lower 
radius of the annulus in the direction of the rejected initial step is begun until a suitable initial iterate 
is found. 


Once the iteration process has successfully started, circles around each iterate are constructed (using 
two theoretical bounds) that are known to bound the root closest to that iterate; the Laguerre step size 
is constantly tested against the radii of these circles and modification of the step is made when it is 
deemed to be too large or when the polynomial value does not decrease in the direction of the step. For 
this reason, the roots are normally found in order of increasing magnitude, thus minimizing the 
roundoff errors resulting from deflation. 


Evaluation of the polynomial and its derivatives at a real iterate is exactly Horner’s method. 
Evaluation at a complex iterate is a modification of Horner’s method that saves approximately half of 
the multiplications. This modification takes advantage of the fact that the Horner recurrence is sym- 
metric with respect to complex conjugation. 


F FOOT uses a sophisticated technique to determine when an approximation Z;, should be accepted as a 
root. As the polynomial is being evaluated at Z,, a bound for the evaluation roundoff error is also being 
computed. If the polynomial value is less than this bound, Z;, is accepted as a root. Z,, can also be 
accepted as a root if the value of the polynomial is decreasing but the size of the Laguerre step has 
become negligible. Before an approximation Z;,, is used in an evaluation, its imaginary part is set to zero 
if this part is small compared to the step size. This improves performance, since real-number 
evaluations are faster than complex evaluations. If the Laguerre step size has become negligible but the 
polynomial is not decreasing, then the message PEOOT #31 lure is reported and the computation 
stops. This is expected never to occur in practice. 
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As the polynomial is being evaluated, the coefficients of the quotient polynomial (by either a linear or 
quadratic factor corresponding to the Z,) are also computed. When an approximation Z, , is accepted as 
a root, this quotient polynomial becomes the polynomial whose roots are sought, and the process begins 
again. 


Multiple Zeros 


No polynomial rootfinder, including FROOT, can consistently locate zeros of high multiplicity with 
arbitrary accuracy. The general rule-of-thumb for FEOOT is that for multiple or nearly-multiple zeros, 
resolution of the root is approximately 12/K significant digits, where K is the multiplicity of the root. 


Accuracy 


PROOT’s criterion for accuracy is that the coefficients of the polynomial reconstructed from the cal- 
culated roots should closely resemble the original coefficients. 


We will illustrate PF:}007T’s performance with isolated zeros using the 100th degree polynomial 
100 

P(Z) = > 2 
k=0 


Of the 200 real and imaginary components of the calculated roots, about half were found to 12 digit 
accuracy. Of the rest, the error did not exceed a few counts in the 12th digit. 


The polynomial (Z + 1)”° with all 20 roots equal to —1 was solved by PRUICIT to yield the following 
roots. 

(—.997874038627,0) 
(—.934656570635,0) 

(—.947080146258, — 160105886062) 
(—.947080146258,.160105886062) 
(—.678701343788, —6.24034855342E — 2) 
(—.678701343788,6.24034855342E — 2) 
(—.815082852233, —.270565874916) 
(—.815082852233,.270565874916) 
(—.725960092383, —.178602450179) 
(—.725960092383,.178602450179) 
(—.934932478844, —.326980158732) 
(—.934932478844,.326980158732) 
(—1.06905713438, —.337946194292) 

(— 1.06905713438,.337946194292) 
(—1.19977533452, —.295162714497) 
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(—1.19977533452,.295162714497) 
(—1.30383056467, —.200016185042) 
(—1.30383056467,.200016185042) 
(—1.3593147483,7.00833934259B — 2) 
(—1.3593147483, — 7.00833934259K — 2) 


The roots appear inherently inaccurate due to the high multiplicity of —1 as a root. Between 0 and 1 
correct digits were expected, even though the first zero found was better than this. However, the re- 
constructed coefficients are very close and are shown below (rounded to 12 digits). 





Original Reconstructed 
Coefficients Coefficients 

1 1 

20 20 

190 190.000000001 

1140 1140 

4845 4845.00000003 

15504 15504 

38760 38760.0000003 

77520 77520.0000007 


125970 125970.000001 
167960 167960.000002 
184756 184756.000002 
167960 167960.000003 
125970 125970.000002 


77520 77520.0000015 
38760 38760.0000009 
15504 15504.0000004 
4845 4845.00000011 
1140 1140.00000004 
190 190.000000042 
20 20.0000000344 


1 1.00000001018 
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Time Performance 


The speed of the FF:O0T function is illustrated in the following table. The times given are those re- 
quired to calculate all the roots of the polynomial 


N 
P(Z) = > 2 
k=0 
for values of N given in the Degree column. 


Note that times are approximate. 





Degree | Time (sec) 














Section 13 


Finite Fourier Transform 





Keyword 


The finite Fourier transform is a key step in solving many problems in mathematics, physics, and 
engineering, such as problems in signal processing and differential equations. 


Given a set of N complex data points Zp, Z,,..., Zy— , the finite Fourier transform will return another 
set of N complex values Wo, W,,..., Wy—,, such that for k = 0, 1,..., N — 1, 


N-1 . : 
Z- > W, (cos 22M + isin aah) 
j=0 N N 


The W’s then represent the complex amplitudes of the various frequency components of the signal 
represented by the data points. The values for the W’s are given by the formula 


— ~2nkj =2nkj 
W,=1/N > Z, (cos 2a + isin aun) 
5 k=0 N N 


This formula holds for any number of data points. The Math Pac uses the Cooley-Tukey algorithm and 
the internal language of the HP-71 to achieve excellent speed and accuracy in the calculation of the 
finite Fourier transform. This requires, however, that N be an integral power of 2; for example, 2, 4, 8, 
16, 32, 64, and 128 are all acceptable values for the number of complex data points. 


To use the finite Fourier transform, store your complex data points Zo,..., Zj—1 as successive elements 
of an N-element complex array with Zy as the first element, Z, as the second element, and so on. Aside 
from the total number of elements in the array, which indicates to the Math Pac the number of com- 
plex data points, the dimensions of the array are irrelevant. For example, each of the following eight- 
element arrays 
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(1,2) 
(3,4) 
(5,6) 
(7,8) 
(9,10) 
(11,12) 
(13,14) 
(15,16) 


(1,2) (3,4) 
(5,6) (7,8) 
(9,10) (11,12) 
(13,14) (15,16) 


(1,2) (3,4) (5,6) (7,8) 
(9,10) (11,12) (13,14) (15,16) 


[«.2) (3,4) (5,6) (7,8) (9,10) (11,12) (13,14) (25,16) | 


can represent the set of input data points 
{(1,2),(3,4),(5,6),(7,8),(9,10),(11,12),(13,14),(15,16) } 


The array in which you wish the transformed data to be stored must be complex type. If the number of 
input data points is N, and if the result array is a vector, it will be redimensioned to have N elements. 
If the result array is a matrix, it will be redimensioned to have N rows and one column. The results of 
the finite Fourier transform Wo,...,W,_; will be returned with the complex values stored in successive 
elements of this N-element complex result array—the same form as the data points. 


The number of data points you can use is limited only by the amount of available memory and by the 
requirement that the number of data points be a non-negative integral power of 2. 
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FOUR Finite Fourier Transform 





MAT W=FOURCZ> 


where Z is an N-element complex array, either a vector or matrix, N is the number of complex data 
points, which must be a non-negative integer power of 2, and W is a complex array, either a vector or 
matrix. 





If W is a vector, implicitly redimensions W to have N elements; if W is a matrix, implicitly redimensions W 
to have N rows and one column. W will be assigned the complex values of the finite Fourier transform of 
the data points represented by Z. 


To halt operation, press twice. 
Not usable in CALC mode. 





Example 


The following example computes the finite Fourier transform of the input data set ((1,2), (3,4), (5,6), 
(7,8), (9,10), (11,12), (13,14), (15,16). 


10 OPTION BASE 1 


20 COMPLEX SHORT A(8),B(1,2) A contains the data set, and B, after 
redimensioning, contains the transform of the 
data. 


30 MAT INPUT A 
40 MAT B=FOUR(A) 
50 MAT DISP B 


RUN 
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3.4142,1.41429 


2,83 
1.4142,-,585 
£219 
38579,-1.41422 





61.4142,-3.41429 





Additional Information 


Time Performance 


The approximate time required by FUUF to return the transform, based on the number of data points, 
is shown in this table. 








Number of | Transform Time 
Data Points (Seconds) 
1 0.07 
2 0.11 
4 0.26 
8 0.75 
16 1.9 
32 47 
64 1 
128 25 
256 55) 
512 120 
1024 260 
2048 558 














Section 13: Finite Fourier Transform 137 


Relation Between the Finite and Continuous Fourier Transform 


The finite Fourier transform is most often used as an approximation to the continuous (infinite) Fou- 
rier transform. To understand in what sense it is an approximation, and to understand the effects of 
various choices to be made in using this approximation, it is most useful to have the direct relationship 
between the continuous and finite transforms. 


If Z(x) is a complex valued function, its continuous Fourier transform is defined to be 
ca . 
Wi) = [~ Z(x) exp(—2nifx) dx 
00 


If we have a set of N complex data points Zp, Z,,..., Zy—; given by sampling the function Z at N 
equally spaced points 


Zp = 2(xo + RAx) fork = 0,1,...,N — 1, 


and then find the finite Fourier transform Wo, W,,..., Wy— of this data set, we can relate these values 
to the values of the continuous Fourier transform W(f) as follows. For k = 0,..., N— 1, 


W, = (r/N) Wk/NAx) where r = exp(—2zixy). 


W is a “smeared” version of the true continuous Fourier transform W. To get W from W, you must 
average W in two important but very different ways. The first type of averaging that occurs can be 
described by defining a new function A(f) intermediate between W and W. 


Af = > Wf + k/Ax) 


k=—o@ 


This says that the value of A at a point f is equal to the sum of the values of W at all points that are 
integer multiples of the limiting frequency 1/Ax away from f. In particular, if W consists of a small 
bump centered at the origin, then A will consist of an infinite sequence of bumps spaced 1/Ax units 
apart. This is the aspect of the finite Fourier transform that gives rise to aliasing: any frequency that 
occurs in W (that is, W has a nonzero value there) will give rise to a nonzero value for A (and also W) 
somewhere in the interval [0, 1/Ax] no matter what the original frequency was. For this reason, you 
should choose Ax small enough so that 1/Ax is larger than the distance between the largest and small- 
est f’s that you suspect will occur in W. Since most functions occuring in actual situations (and all 
real-valued functions) have continuous Fourier transforms that are roughly symmetric about the origin, 
if a frequency fg occurs in W, it is likely that —f also occurs in W. For the finite Fourier transform to 
contain both frequencies without aliasing, 1/Ax must be larger than 2f/). If we define the largest fre- 
quency occuring in W as Af, we can express the no-aliasing requirement as AfAx < 1/2. 
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The second type of averaging that occurs when going between W and W is much more local in nature 
than the first. It results in a loss of frequency resolution in W as compared with W; more precisely, 


Wi) = (Nax) [™ sine(gNAx) A(f — 2) dg 


lifa=0, 
where sinc(a) = 5 
Sin (ra) otherwise. 
7a 


Since sinc (gNAx) consists primarily of a bump with width inversely proportional to NAx, W is more 
blurred (compared to W) for smaller values of NAx. This is not a serious problem unless W has a large 
value at a frequency that is not a multiple of the fundamental frequency N/Ax. In this case, the “side 
lobes” of the sinc function become evident in W. This can be reduced somewhat by multiplying the 
data values Z;, by a smoothing function G(k) before taking the finite Fourier transform. This results in 
an averaging function that has smaller side lobes than the sinc function. One example of such a func- 
tion is the Hanning function G(k) = (1/2)(1 — cos(2rk/N)). 


Inverse Finite Fourier Transform 


Many applications of the finite Fourier transform involve taking the transform of a set of data points, 
operating on the transformed values (for example, increasing or decreasing the amplitudes), and then 
retransforming the data using the inverse Fourier transform defined by 


2, = "5, W; (cos 228 + i sin 2281) 
i — es SD 2 pnp 


You can also use the F (LIF! keyword to compute the inverse finite Fourier transform in a simple way. If 
W is an N-element complex array for which you want the inverse finite Fourier transform: 


1. Redimension W to have N rows and one column (if W is an array with only one column, then no 
redimensioning is necessary). 


2. Take the transpose (T FH) of W. This produces the complex conjugate of W, without changing the 
order of the elements. 


3. Take the finite Fourier transform of the result. 


4. Take the transpose of the result of the finite Fourier transform and scalar multiply this result by 
N. This will produce the inverse finite Fourier transform of the original array. 
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Example. This illustrates an application of the finite Fourier transform, and shows the procedure for 
obtaining the inverse finite Fourier transform. 


Suppose we want to find the steady state solution Z(x) of the inhomogeneous differential equation 
Z"(x) + 3Z(x) + 12Z(x) = P(x) 


where P(x) is a function for which we have sampling data. If we denote the (continuous) Fourier trans- 
form of any function Q by Q, by taking the Fourier transform of the above equation we arrive at 


—f?2(f) + 3ifZ( + 122() = Pip). 


Solving this equation algebraically we obtain 


yet PGs 
20 (—f? + 12) + 3if 


If we can get a good approximation of P, we can easily calculate the right hand side of this equation. 


From this result we can obtain the solution to the original equation by taking the inverse Fourier 
transform. 


For simplicity, we will assume that the equation has been scaled so the P(x) has unit period, and that 
the highest frequency component of P is (approximately) 30 times the fundamental frequency. Sam- 
pling P 64 times in one period will then suffice to avoid aliasing. 


Rather than prompt the user for 64 complex data points representing the sampling of P, the program 
below uses a relatively simple function for P, although you could use values from any other source 
equally well. 


10 OPTION BASE 1 

20 COMPLEX P(64),Q(64,1),Z(1,64) P will contain the data points representing the 
sampling of P. Q will represent P and eventually 
P/(—f? + 3if + 12). Z will represent the solution 
to the differential equation. 

30 COMPLEX T T is a complex scalar for use in the complex 
division. 

40 DISP “Working; please wait.” 

50 RADIANS 
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60 FOR I=1 TO 64 

70 R=PI*1/32 

80 P(l)=( 6000*COS(3*R)*SIN(7.5*R)* 
COS(5.5*R) , 4000*COS(13*R)+ 
3500*SIN(11%*R) ) 


90 NEXT | 
100 MAT Q=FOUR(P) 
110 FOR F=—31 TO 32 


120 J=MOD(F,64)+1 


130 T=(—F°2+12,3%*F) 


140 Q(J,1)=Q(J,1)/T 
150 NEXT F 
160 MAT Q=TRN(Q) 


170 MAT Z=FOUR(Q) 
180 MAT Z=TRN(Z) 


190 MAT Z=(64)*Z 
200 COMPLEX Z(64,1) 
210 DISP “The values are” 


220 MAT DISP USinwG 
“X,C(MDDD.D,MDDD.D)";Z 


This is the sampling routine that assigns to P 
the values of the complex-valued functions 
represented by the right-hand side of line 80, 
sampled at 64 equally spaced points. 


Q now represents P. 


F represents the frequency variable and spans 
the full range of frequencies, positive and 
negative, that we expect to occur in P. 


J represents the number of the element in the Q 
array where the amplitude of the frequency F is 
stored. 

T will be the denominator of the complex 
fraction. 


Q now represents P/(—f? + 3if + 12). 

This starts the procedure that assigns the values 
of the inverse Fourier transform to Z. The 
transpose is used here to take the conjugate of Q. 


The transpose is used here for conjugation as 
well. 


The values displayed will represent the complex 
values of the steady state solution of the 
differential equation sampled at 64 equally 
spaced points in one period. 
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Fourier Sine/Cosine Series 


There is another transform closely related to the finite Fourier transform that is applicable when the 
data points Z;, are purely real (that is, their imaginary parts are equal to zero). This is the Fourier 
series transformation, which takes a set of 2N (real) data points Zp, Z;, ..., Z2y—, and returns a set of 
2N + 1 real values Ap, Aj, ...,Ay, By, -.., By with the property that 


Ao Ny 2ajk « Qnjk 
Z, = 2 + 2 A; 000 oy + B; sin “I 
If Wo, Wj, -.-, Won—1 are the complex values of the fir te Fourier transform of the real data points 


Zo, ---, Z,n—1, then the Fourier series values are given by 
A; = 2REPT(W;) for j = 0, ..., N—1, 


Ay = REPT(Wy) 


By = —2IMPT(W,) for j = 1,..., N. 





Appendix A 
Owner’s Information 


Installing and Removing the Math Pac Module 


The math module can be plugged into any one of the four ROM ports on the front edge of the 
computer. 





CAUTIONS 


e Be sure to turn off the HP-71 (press [f] [OFF]) before installing or removing the module. 


« If you have removed a module to make a port available for the math module, before installing the 
math module, turn the computer on and then off to reset internal pointers. 


« Do not place fingers, tools, or other objects into any of the ports. Such actions could result in minor 
electircal shock hazard and interference with pacemaker devices worn by some persons. Damage 
to port contacts and internal circuitry could also result. 


e If a module jams when inserted into a port, it may be upside down. Attempting to force it further 
may result in damage to the computer or the module. 

e Handle the plug-in modules very carefully while they are out of the computer. Do not insert any 
objects in the module connecter socket. Always keep a blank module in the computer port when a 
module is not installed. Failure to observe these cautions may result in damage to the module or 
the computer. 








Limited One-Year Warranty 


What We Will Do 


The Math Pac is warranted by Hewlett-Packard against defects in materials and workmanship affect- 
ing electronic and mechanical performance, but not software content, for one year from the date of 
original purchase. If you sell your unit or give it as a gift, the warranty is transferred to the new owner 
and remains in effect for the original one-year period. During the warranty period, we will repair or, at 
our option, replace at no charge a product that proves to be defective, provided you return the product, 
shipping prepaid, to a Hewlett-Packard service center. 
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What Is Not Covered 


This warranty does not apply if the product has been damaged by accident or misuse or as the result of 
service or modification by other than an authorized Hewlett-Packard service center. 


No other express warranty is given. The repair or replacement of a product is your exclusive remedy, 


ANY OTHER IMPLIED WARRANTY OF MERCHANTABILITY OR FITNESS IS LIMITED 
TO THE ONE-YEAR DURATION OF THIS WRITTEN WARRANTY. Some states, provinces, 
or countries don’t allow limitations on how long an implied warranty lasts, so the above limitation may 
not apply to you. IN NO EVENT SHALL HEWLETT-PACKARD COMPANY BE LIABLE 
FOR CONSEQUENTIAL DAMAGES. Some states, provinces, or countries do not allow the exclu- 
sion or limitation of incidental or consequential damages, so the above limitation may not apply to you. 


This warranty gives you specific legal rights, and you may also have other rights which may vary from 
state to state, province to province, or country to country. 


Warranty for Consumer Transactions in the United Kingdom 


This warranty shall not apply to consumer transactions and shall not affect the statutory rights of a 
consumer. In relation to such transactions, the rights and obligations of Seller and Buyer shall be 
determined by statute. 


Obligation To Make Changes 


Products are sold on the basis of specifications applicable at the time of manufacture. Hewlett-Packard 
shall have no obligation to modify or update products once sold. 


Warranty Information 


If you have any questions concerning this warranty, please contact an authorized Hewlett-Packard 
dealer or a Hewlett-Packard sales and service office. Should you be unable to contact them, please 
contact: 


e In the United States: 


Hewlett-Packard Company 
Personal Computer Group 
Customer Communications 
11000 Wolfe Road 
Cupertino, CA 95014 


Toll-Free Number: (800) FOR-HPPC (800 367-4772) 
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e In Europe: 
Hewlett-Packard S.A. 
150, route du Nant-d’Avril 
P.O. Box CH-1217 Meyrin 2 
Geneva 
Switzerland 
Telephone: (022) 83 81 11 


Note: Do not send products to this address for repair. 


e In other countries: 


Hewlett-Packard Intercontinental 
3495 Deer Creek Rd. 
Palo Alto, CA 94304 
U.S.A. 
Telephone: (415) 857-1501 


Note: Do not send products to this address for repair. 


Service 


Service Centers 


Hewlett-Packard maintains service centers in most major countries throughout the world. You may 
have your product repaired at a Hewlett-Packard service center any time it needs service, whether the 
unit is under warranty or not. There is a charge for repairs after the one-year warranty period. 


Hewlett-Packard computer products normally are repaired and reshipped within five (5) working days 
of receipt at any service center. This is an average time and could vary depending on the time of year 
and work load at the service center. The total time you are without you product will depend largely on 
the shipping time. 
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Obtaining Repair Service in the United States 


The Hewlett-Packard United States Service Center for battery-powered computational devices is lo- 


cated in Corvallis, Oregon: 


Hewlett-Packard Company 
Service Department 
P.O. Box 999 


Corvallis, OR 97339, U.S.A. 


or 


1030 N.E. Circle Blvd. 


Corvallis, OR 97330, U.S.A. 


Telephone: (503) 757-2000 


Obtaining Repair Service in Europe 


Service centers are maintained at the following locations. For countries not listed, contact the dealer 


where you purchased your unit. 


AUSTRIA 

HEWLETT-PACKARD Ges.m.b.H. 
Kleinrechner-Service 
Wagramerstrasse-Lieblgasse 1 
A-1220 Wien (Vienna) 

Telephone: (0222) 23 65 11 


BELGIUM 

HEWLETT-PACKARD BELGIUM SA/NV 
Woluwedal 100 

B-1200 Brussels 

Telephone: (02) 762 32 00 


DENMARK 
HEWLETT-PACKARD A/S 
Datavej 52 

DK-3460 Birkerod (Copenhagen) 
Telephone: (02) 81 66 40 


EASTERN EUROPE 
Refer to the address listed under Austria. 


FINLAND 
HEWLETT-PACKARD OY 
Revontulentie 7 

SF-02100 Espoo 10 (Helsinki) 
Telephone: (90) 455 02 11 


FRANCE 

HEWLETT-PACKARD FRANCE 
Division Informatique Personnelle 
S.A.V, Calculateurs de Poche 
F-91947 Les Ulis Cedex 
Telephone: (6) 907 78 25 


GERMANY 
HEWLETT-PACKARD GmbH 
Kleinrechner-Service 
Vertriebszentrale 

Berner Strasse 117 
Postfach 560 140 

D-6000 Frankfurt 56 
Telephone: (611) 50041 


ITALY 

HEWLETT-PACKARD ITALIANA S.P.A. 
Casella postale 3645 (Milano) 

Via G. Di Vittorio, 9 

|-20063 Cernusco Sul Naviglio (Milan) 
Telephone: (2) 90 36 91 


NETHERLANDS 

HEWLETT-PACKARD NEDERLAND B.V. 
Van Heuven Goedhartlaan 121 

N-1181 KK Amstelveen (Amsterdam) 
P.O. Box 667 

Telephone: (020) 472021 


NORWAY 

HEWLETT-PACKARD NORGE A/S 
P.O. Box 34 

Oesterndalen 18 

N-1345 Oesteraas (Oslo) 
Telephone: (2) 17 11 80 


SPAIN 

HEWLETT-PACKARD ESPANOLA S.A. 
Calle Jerez 3 

E-Madrid 16 

Telephone: (1) 458 2600 


SWEDEN 

HEWLETT-PACKARD SVERIGE AB 
Skalholtsgatan 9, Kista 

Box 19 

S-163 93 Spanga (Stockholm) 
Telephone: (08) 750 20 00 


SWITZERLAND 
HEWLETT-PACKARD (SCHWEIZ) AG 
Kleinrechner-Service 

Allmend 2 

CH-8967 Widen 

Telephone: (057) 31 21 11 


UNITED KINGDOM 
HEWLETT-PACKARD Ltd 
King Street Lane 
GB-Winnersh, Wokingham 
Berkshire RG11 5AR 
Telephone: (0734) 784 774 
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International Service Information 


Not all Hewlett-Packard service centers offer service for all models of HP products. However, if you 
bought your product from an authorized Hewlett-Packard dealer, you can be sure that service is avail- 
able in the country where you bought it. 


If you happen to be outside of the country where you bought your unit, you can contact the local 
Hewlett-Packard service center to see if service is available for it. If service is unavailable, please ship 
the unit to the address listed above under Obtaining Repair Service in the United States. A list of 
service centers for other countries can be obtained by writing to that address. 


All shipping, reimportation arrangements, and customs costs are your responsibility. 


Service Repair Charge 


There is a standard repair charge for out-of-warranty repairs. The repair charges include all labor and 
materials. In the United States, the full charge is subject to the customer's local sales tax. 


Computer products damaged by accident or misuse are not covered by the fixed repair charge. In these 
cases, repair charges will be individually determined based on time and materials. 

Service Warranty 

Any out-of-warranty repairs are warranted against defects in materials and workmanship for a period 
of 90 days from date of service. 

Shipping Instructions 


Should your product require service, return it with the following items: 
* A completed Service Card, including a description of the problem. 


° A sales receipt or other documentary proof of purchase date if the one-year warranty has not 
expired. 


The product, the Service Card, a brief description of the problem, and (if required) the proof of pur- 
chase date should be packaged in adequate protective packaging to prevent in-transit damage. Such 
damage is not covered by the one-year limited warranty; Hewlett-Packard suggests that you insure the 
shipment to the service center. The packaged product should be shipped to the nearest Hewlett- 
Packard designated collection point or service center. Contact your dealer for assistance. 
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Whether the product is under warranty or not, it is your responsibility to pay shipping charges for 
delivery to the Hewlett-Packard service center. 


After warranty repairs are completed, the service center returns the product with postage prepaid. On 
out-of-warranty repairs in the United States and some other countries, the product is returned C.O.D. 
(covering shipping costs and the service charge). 


Further Information 


Service contracts are not available. Computer products circuitry and design are proprietary to Hewlett- 
Packard, and service manuals are not available to customers. Should other problems or questions arise 
regarding repairs, please call your nearest Hewlett-Packard service center 


When You Need Help 


Hewlett-Packard is committed to providing after-sale support to its customers. To this end, our cus- 
tomer support department has established phone numbers that you can call if you have questions about 
this product. 


Product Information. For information about Hewlett-Packard dealers, products, and prices, call the 
toll-free number below: 


(800) FOR-HPPC 
(800 367-4772) 


Technical Assistance. For technical assistance with your product, call the number below: 
(503) 754-6666 
For either product information or technical assistance, you can also write to: 


Hewlett-Packard Company 
Personal Computer Group 
Customer Communications 
11000 Wolfe Road 
Cupertino, CA 95014 





Appendix B 
Memory Requirements 


The Math Pac reserves 43.5 bytes of read/write memory for its own uses. In addition, small amounts of 
memory are temporarily used for routine overhead purposes. Significant amounts of memory can be 
used to declare complex variables and arrays (see page 20), and to redimension arrays to a larger size, 
but this memory usage is easily determined. This appendix lists the amounts of tempory memory used 
by other Math Pac operations. 





Item Memory Required For Operation 


Matrix operations 
OETSA? 2N(4N + 1) bytes, where A is an N x N matrix. 


MAT PRINT USING | 14 bytes. 
MAT DISP USING 14 bytes. 


MAT IHPUT 40 bytes. 

MAT A=AtA Requires additional memory only if an operand array is used for the result 
MAT A=AtB array. If the product (that is, the redimensioned array A) is M x N (for vec- 
MAT A=BHA tors, let N = 1), then the memory required is: 


3MN bytes, if A is type INTEGER. 

4.5MN bytes, if A is type SHORT. 

8MN bytes, if A is type FEAL. 

QMN bytes, if Ais type COMPLEX SHORT, 
16MN bytes, if A is type COMPLEX. 


NAT A=TRHC ADEA Requires additional memory only if an operand array is used for the result 
MAT A=TRHCAD# array. If the product (that is, the redimensioned array A) is M x N (for vec- 
MAT A=TRH( BDA tors, let N = 1), then the memory required is: 

3MN bytes, if A is type INTEGER. 

4.5MN bytes, if Ais type SHORT. 

8MN bytes, if A is type REAL. 

QMN bytes, if Ais type COMPLEX SHORT. 

16MN bytes, if A is type COMPLEX. 
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MAT A=TRN¢A> 


MAT B=PROOTCA? 


MAT B=FOURCA> 


FHROOT 


INTEGRAL 








Item te Memory Required For Operation 

MAT B=IHVcAd AisN XN. 
If Ais REAL, SHORT or INTEGER and Bis FEAL: 

4N bytes. 


If Ais REAL, SHORT or INTEGER and B is not REAL: 
4N(2N + 1) bytes. 

If Ais COMPLEX or COMPLEX 
8N(4N + 1) bytes. 


Ais N x N and B is N x P (for vectors, let P = 1). 
If Ais REAL, SHORT, or INTEGER and B is REAL, SHORT, or 
INTEGER: 
4N(2N + 4P + 1) bytes. 
lf A is REAL, SHORT, or INTEGER and B is COMPLEX or 
COMPLEX SHORT: 
4N(2N + 8P + 1) bytes. 
If Ais COMPLEX or COMPLEX SHORT: 
8N(4N + 4P + 1) bytes. 
If Ais M x Nand INTEGER: 
MN/2 bytes. 
If operand and result matrix are different, or if A is not IHTEGER, no extra 
memory is needed. 
A represents an Nth degree polynomial. 
21N + 261 bytes. 


A contains N elements. 
If Bis COMPLEX SHORT: 
16N bytes. 
If Bis COMPLE type, no extra memory is needed. 
112.5 bytes if FHROOT is not nested. 
96.5 additional bytes for each level of nesting. 
208.5 bytes if INTEGRAL is not nested. 
192.5 additional bytes for each level of nesting. 





SHORT: 
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Error Conditions 


The Math Pac reports two classes of error messages. 


° Math Pac error messages. These have a LEX ID number of 2. These error messages are explained 
in the first table. 


° HP-71 error messages that are reported by Math Pac functions. These have a LEX ID number of 
0. These error messages are explained in the second table. 


Math Pac Error Messages 





pNumber 


1 








Error Message and Condition 


#01MS 
© DOTA, B>: A or B is a matrix. 


° DETCA?, MAT B=INV¢CAD, MAT B=TRN¢CA?, MAT A=I0h, 
MAT X=S'S¢A,Y3: A or B is a vector. 


°© MATA=I0H¢/2: only one redimensioning subscript specified. 
© MAT A=operation « operand array(s)\: number of subscripts of A is not the same as the 
number of subscripts required for the result of the operation. 
Hot Square 
* DET¢A®, MAT A=10H, MAT B=INW CAD, MAT X=S'S¢A,B3: Ais a matrix but 
the number of rows of A is not equal to the number of columns. 
© MAT A=IONC,f2:1 ¥ j. 








151 





a 
P<} 


Append 








Number 


Error Message and Condition 








Conformability 
e MAT A=B+C, MAT A=B-C: B and C are not conformable for addition (the number of 
rows are unequal or the number of columns are unequal) 


e MAT A=B#C: B and C are not cdnformable for multiplication (B is a vector or the num- 
ber of columns of B is not equal to the number of rows of C). 

e MAT A=TEH¢Bs*C: B and C are not conformable for transpose multiplication (B is a 
vector or the number of rows of B is not equal to the number of rows of C). 


e MAT A,B»: Although A is a square matrix, A and B are not conformable for 
multiplication. 





e OOT«A, B®: Although A and B are vectors, the number of elements of A is not equal to 
the number of elements of B. 
Parameter Redim 
e The result array of a HT statement is a subprogram parameter. The statement requires 
array redimensioning, which changes the number of array elements. 
Hesting Error 
e More than five FHROOT or IHTEGRAL keywords are nested. 
Kybd FA in FAROOTs INTEGRAL 
e Attempting to execute FHROOT or INTEGRAL from the keyboard in BASIC mode, 
and the function whose root or integral is sought is a user-defined function. 


e Attempting to execute a user-defined function from the keyboard while an FHF: CT or 
IHTEGRAL execution is suspended during the evaluation of the function whose root or 
integral is sought. 


Function Interrupted 


e Interrupting DET* A>, CHORMN¢ A), RMORMCA?, FHORMCA®, or COTA, B® by 
pressing [ATTN] twice. 














Bad Array Size 
e MAT B=FOUR«A® where the number of elements of A is not a non-negative integral 
power of two. 


e MAT B=PROOT «A> where A has only one element. 
PROOT Failure 
e FROOT failed to find a root. 
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Number Error Message and Condition 





10 GANMA=InF 

e GAMMACX> where X is a non-positive integer. 
1 ATFHHE +-13 

© ATAWHE 19 of ATANHE-13 


No Initislizatioan 
ee © The Math ROM cannot initialize due to insufficient memory. This ROM requires 43.5 


bytes of user memory for its own use. This memory must be available before plugging in 
the module. 









HP-71 Error Messages 








s 
Number Error Message and Condition 
11 Invalid Arg 
© BVAL“B$,R>, BSTR#¢X,R: The rounded integer value of R is not equal to 2, 8, or 
16. 


° BYAL¢B$,R>: B$ is not a valid string representation of a number in base A. 

© BSTR#«X,A: The rounded integer value of X is not in the interval [0,1E12). 

° BVAL“B$.R>: The decimal equivalent of B$ exceeds 999,999,999,999. 

e LBMO*CA,N2, UBHO¢A,N>: The rounded integer value of N is not equal to 1 or 2. 

° An illegal subscript has been used ina MAT COM, MAT IOGH, MAT ZER, 
COMPLE, or COMPLEX SHORT statement. 

24 Insufficient Memory 

© Appendix B gives the memory requirements for various Math Pac operations. 

31 Data Type 

© A scalar (real or complex) has been used where an array is required or vice-versa. 


© A complex ‘ype (scalar or array) has been used where a real type (scalar or array) is 
required or vice-versa. 

















| Number 


Error Message and Condition 


32 Ho Data 


e Attempting to execute DETL before the first completion of MAT... IN! with a real- 
type argument or MAT... S with a real-type first argument. 
e Attempting to execute F'VALUE or FGU 
FHROOT keyword. 
e Attempting to execute IVVALUE or IBOUHD before any INTEGRAL keyword has 
completed the first evaluation of the function whose integral is sought. 
» Attempting to execute F'/AR while no FHROOT is evaluating the function whose root is 
| sought. 


» Attempting to execute I'\'AF while no INTEGRAL is evaluating the function whose 
integral is sought. 








before the first completion of an 


46 Invalid USING 
e Formatting a real expression with a complex IMAGE field or vice-versa. 
79 Tllegsl Context 


e Attempting to execute INTEGRAL or FHROOT from CALC mode in any way except 
by direct execution. 


80 Invalid Parameter 


« MAT INPUT attempts to execute an expression in the MAT INFLUT response line 
where that expression calls a user-defined function. 
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Attention Key Actions 


The way operates during the execution of each the following three keywords is described on the 
referenced page. 


MAT INPUT Refer to page 54. 
FHROOT Refer to page 97. 
INTEGRAL Refer to page 111. 


The keywords listed below in this appendix can be aborted by pressing the key once or twice. 


Array Output Statements 


All Math Pac array output statements (MAT OCISF/PRIHT[USIMG]) can be halted at any time by 
pressing once. 





Other |: 


The following iit statements may be halted at any time by pressing twice. 


| Statements 


MAT result = [-] operand 


MAT result = operand +/—/# operand 


MAT result = « scalar > [+ operand] 


i 


MAT result IMM! ¢ operand * 





MAT result = 3 operand , operand » 
NAT result = TRH« operand 2[%# operand ] 
MAT result = FOURS® operand » 


i 


MAT result FROOT* operand > 
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Suppose a lengthy program contained a MAT IH‘! statement. Suppose further that you wished to 
abort this program. You press once, and the program does not halt (the SUSP annunciator does 
not turn on). This tells you that the MAT IH‘! statement may be executing, and gives you a chance to 
wait for the result of this MAT IHN! execution, or to abort the MAT IN! execution and the program 
immediately by pressing a second time. In this way the “press N] twice” rule gives a user 
more control over program and statement suspension. 






Pressing once during execution of MAT IH! would suspend the program in the usual way after 
this statement is completed. 


Scalar-Valued Array Functions 





The following scalar-valued array functions can be halted at any time by pressing [ATTN] twice. 











GET« operand > 

OOT« operand , operand >} 
FHORIM< operand > 
CHORM © operand 
EHORM« operand > 


The benefits provided by this “press twice” rule are the same as those described above. However, 
only an error can halt the execution of an expression, so when you press [ATTN] twice to halt any of the 
above functions, the HP-71 will display the error message Furnnction Interrupted. 
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Numeric Exceptions and the IEEE Proposal 


Introduction 


This appendix will discuss IEEE exception handling by Math Pac functions and operations, including 
computation with H3M and Inf arguments, exception flag setting, handling of out-of-range arguments, 
error or warning messages, and default values for I\’L and [i\'2 exceptions. The HP-71 reference 
manual discusses the IEEE proposal for handling math exceptions. Math Pac functions, when appro- 
priate, will set the exception flags IL, OZ, OF, UNF, and INX and report errors or warnings (with 
default results returned) according to the TAF settings for each of these flags. You can refer to the 
appropriate sections of this manual for definitions and/or computational formulas for many of the 
functions described here. 


No exception flags are set by any of the keywords in sections 2 or 3 of this manual, or by Math Pac 
keywords TYPE, ~ (negation of complex numbers), COHJ, COM, IOQHN, ZER, 
MAT OISF/PRINT[USING), LEND, UBND, DETL, FYAR, FYALUE, FGUESS, IVAR, IVALUE, 
and IBOUNO. Remember that exception flags IM, OF, and UNF may be set when values are 
rounded to fit the destination type, such as, for example, assigning {MAMREAL,MAXREAL® toa 
COMPLE SHORT variable or executing HT A=B where A is INTEGER type and B contains ele- 
ments greater than 99999. 





Aside from exceptions occurring during rounding, the statements MAT A=B, MAT A=~B, 
MAT A=TRN<B>, and MAT A=<X> set only the I‘’L exception flag (reporting message 
Signaled Op) and only when Ais INTEGER type and either B contains, or ~ is, a signaling Hal. 
This is because I!/TEGER variables can contain only quiet, not signaling, !/ ats. The same applies to 
MAT IHPUT. 


The cases given for each of the keywords in the tables which follow are evalucted in order from top to 
bottom. 


Note: Throughout this appendix, * represenis any argument. 
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Real Scalar Functions 


These functions are described in section 4 of this manual. Any signaling 4! 3! argument sets I\!L and 
reports message Signaled Op; if TRAP <IL> = 2, then this NaH becomes quiet and the opera- 
tion can continue. With the exception of the 43H function, any quiet 4a} argument returns a Nal! 
result with no exception flags set. (Aside from signaling 12} arguments, the functions IF. 0/JHO and 
4aH# set no exception flags). 


Real Hyperbolic Sine (S1HH*X") 








Argument X Result 
+Inf X; no exception flags set. 
+0 X; no exception flags set. 
* IH® set; UNF, OF set as appropriate. 





Real Hyperbolic Cosine (C05H«X") 








Argument X Result 
+Inf |X|; no exception flags set. 
+0 1; no exception flags set. 


* IHS set; O'/F set as appropriate. 


ie : = 


Real Hyperbolic Tangent (TANH*X») 








Argument X Result 
+Inf SGH¢X3; no exception flags set. 
+0 x 
* IH® set; HF set as appropriate. 





Real Hyperbolic Arc Sine (AS IHH«X>) 





Argument X Result 





+Inf X; no exception flags set. 
+0 X; no exception flags set. 
IH set; LMF set as appropriate. 
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Real Hyperbolic Arc Cosine (ACO SH«X*) 








Argument X Result 
Inf X; no exception flags set. 
x<1 IL set; NH result; message Invalid Arg. 
1 0; no exception flags set. 
* IN# set. 





Real Hyperbolic Arc Tangent (A TAHH«X*) 





Argument X Result 





1 IL set; NaN result; message Invalid Arg. 
1 DW2 set; message ATANH¢ +-13. 
SGM&X> x Inf result if TRAP COY 
SGN¢X) x MAXREAL result with IHX 






mg, 
set if TRAPCOWZ> = 1. 











+0 X; no exception flags set. 
* IN set; UHF set as appropriate. 
Base 2 Logarithm (LOG2¢X) 
Argument X Result 
Inf X; no exception flags set. 
x<0 IVL set; NaH result; message LOG¢neg>. 
+0 OWS set; message LOG (G3, 


—Inf result if TRAPCOVZ> = 2. 

—MAXREAL result with IHX set if TRAPCOVZ> = 1. 
1 0; no exception flags set. 

* IH set. 




















Gamma Function (GAMNMA®X2) 
Argument X Result 
Inf X; no exception flags set. 
+0 OZ set; message GANMA=I1HF. 
CLASS¢X> x Inf result if TRAP COW = 2. 
CLASS¢X> x MAXREAL result with IHS setif TRAP COVES) = 1. 
X < 0 and integral OZ set; message GAMNMA=IHF. 
| —Inf result if TRAPCOVZ> = 2, 
i —MAXREAL result with INX set if TRAPS OVE) = 1, 
* IH set for all X not in the set 











{1, 2, ..., 18}; UHF, OF set as appropriate. 
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Nearest Machine Number (HE IGHEOGR<X, ¥3) 





Arguments 


Result 





X; UNF, IH& set if TRAPSUNF> # 2 andO < |X| < EPS. 
Y; no exception flags set. 
Y; no exception flags set. 








SGN¢X> x MAXREAL; no exception flags set. 

SGNCY) x MINREAL; UHF, INH set if TRAPCUMF > # 2, 
MIMREAL +0 0; no exception flags set. 
—MIHREAL | +0 —0; no exception flags set. 
* * UHF, IHX set if |NEIGHBOR<X, Y3|< EPS and TRAP CUNF> # 2, 





Power of Ten Scaling (SCALE1&°¢X,N*) 

















Arguments 
Result 
N 
: L 
* non-integer | I\/L set; NaH result; message Inwalid Arg. 
+1nf =—Tnft IVL set; HaM result; message Inf #4. 
0 Inf IVL set; HaH result; message Inf ¥ 
eit * X; no exception flags set. 
* —Inf SGM«X) x 0; no exception flags set. 
* Inf SGH¢X3 x Im; no exception flags set. 
* * IH, OF, UNF set as appropriate. 








Complex Functions and Operations 


These functions are described in section 5 of pts manual. For extensions of HP-71 and Math Pac 
functions to complex arguments (4+, -, #, +, *, LOG, EXP, SIM, COS, TAN, SINH, COSH, TANH, 
27, SGN, ABS, 7, and #), er the complex case is discussed here. For the functions 


POLAR, RECT, ARG, and PRO.J, computation at a real argument X is equivalent to computation at 
the complex argument (X, 0). 





Any signaling 4H argument (including real and imaginary parts of complex arguments) sets I \!L. and 
reports message Signaled Op;if TRAP¢IL > = 2, then this H3H becomes quiet and the opera- 
tion can continue. In the following discussion, all references to 13s are to quiet Hats. 
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The following terms are used: 
e Complex denotes complex DATA type. 
e Real denotes real DATA type (e.g., (3, 0) is complex and 3 is real). 
e CaM denotes any complex number with at least one Hah component. 


e CIrif denotes any complex number whose magnitude is If; that is, any complex number with at 
least one +Im# component. 


e CZERO denotes any complex number whose magnitude is 0. 

e Arg(Z) denotes the argument of Z, that is, the infinitely precise value of the Math Pac function 
ARGEZ3. 

« |Z| denotes the magnitude of Z. 

e The complex variables Z and W will also be denoted by (x, y) and (u, v) respectively. 


+, ~ (Addition and Subtraction) 


For real a and complex Z,a + Z = (a + x, y) and Z + a = (x + a, y). For complex Z and W, Z + W 
= (x tu, y + v). IVL is set and message Inf-Inf is reported if any componentwise addition or 
subtraction is equivalent to Inf — Inf; a NaH is returned for the corresponding result component. 
Otherwise, IM, OWF, and LINF are set for each result component as appropriate. 


(Multiplication) 


For real a and complex Z,a x Z = Z x a = (ax, ay). I'L is set and message If #4 is reported if 
any componentwise multiplication is equivalent to (+Inf) x (+0); a Hal is returned for the 
corresponding result component. Otherwise, IM®, 0\'F, and UHF are set for each result component as 
appropriate. 
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For complex Z and W, Z x W is given by the table below. 


Complex x Complex Multiplication (Z - W) 


CHal 
* 
CiInf 
CZERO 
CInf 
* 
* 





Result 
















“Hah, MaN3; no exception flags set. 

“Hah. MaH2; no exception flags set. 

IVL set; “NaH,NaN> result; message Inf #4. 

IVL set; (NaH, NaN result; message Inf #8. 
RECT¢¢Inf, Arg(Z) + Arg(W)> >; no exception flags set. 
RECTCCIn#, Arg(Z) + Arg(W)> >; no exception flags set. 


(xu — yv, xv + yu); IN®, OF, UNF set for each result component as 
appropriate. 





» (Division) 


For real a and complex Z, Z/a = (x/a, y/a). I\‘L is set and message &. is reported if any 
componentwise division is equivalent to (+0)/(+0); a NaH is returned for the corresponding result 
component. IL is set and message Inf .“Inf is reported if any componentwise division is equivalent 
to (+Inf)/(+Imf); a Hah is returned for the corresponding result component. [!'2 is set and mes- 
sage “Zero is reported if any componentwise division is equivalent to 7/(+0) where T' is neither a 
Hal, Inf, or +0; Inf of the appropriate sign is returned for the corresponding result component if 
TRAP ¢OM2» = 2; MAXREAL of the appropriate sign is returned with IH™ set for the corresponding 
result component if TRAP < O23 = 1, Otherwise, IM, OF, and UNF are set for each result compo- 
nent as appropriate. 








For complex Z, we define the following. If Z = CZERO, then 1/Z is defined to be (CLAS 
-SGH(y)). If Z = CInf, then 1/Z is defined to be (SGN(x) x 0, -SGH(y) x 0). 





x) x Inf, 
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For real a and complex Z, a/Z is given by the table below. 


peal Complex Division (a / Z) 


Arguments 

a Zz 
Hal * 
* CHaH 
Lit CiInf 
+0 CZERO 
+Inf | CZERO 
* CZERO 
* CiInf 
Tint * 
* * 


| (a/|Z|*) x COMI¢Z3 (real x complex multiplication); IH 


Result 
©HMat. Nah; no exception flags set. 
“NaN. Mah; no exception flags set. 
IML set; (HaN, NaN? result; message Inf «Inf 
IVL set; (HaN.Nah3 result; message 4.4. 


SGN(a) x (1/Z) (real x complex multiplication); no exception flags set. 





OWZ set; message “Zera. 
a x (1/Z) (real x complex multiplication) result if TRAP «Oe > 

a xX (1/Z) (real x complex multiplication) result with +Imf result component | 

replaced by +MAXREAL and IHx set if TRAP COVWe2 > 


a X (1/Z) (real x complex multiplication); no exception flags set. 

a x COMIEZ> (real x complex multiplication); I'L set and message Inf #1 | 
reported if any componentwise multiplication is equivalent to (+ Inf) x (+0); a 
NaH is returned for the corresponding result component. Otherwise, no excep- | 
tion flags set. 








IWF, LMF set for 





each result component as appropriate. 
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For complex Z and W, W/Z is given by the table below. 


Complex / Complex Division (W/Z) 




















Arguments 
Result 
Ww Zz 
CHaH * “Hal, Ha; no exception flags set. 
*” CHa ©Mah, Mah>; no exception flags set. 
CZERO | CZERO | I'!L set; ¢N3H,NaH> result; message 4.4. 
CInf CInf IVL set; “NaN. NaN? result; message Inf-Inf. 
CiInf CZERO | W x (1/Z) (complex x complex multiplication); no exception flags set. 
* CZERO | 0\2 set; message “Zero. 
W x (1/Z) (complex x complex multiplication) result if TRAP © OW29 = 2. 
W x (1/Z) (complex x complex multiplication) result with + Irif result compo- 
nent(s) replaced by +MAXREAL and IN setif TRAPCOVWEs = 1. 
* CiInf W x (1/Z) (complex x complex multiplication); no exception flags set. 
* * (W x CONJSZ>)/|Z|2? (complex x complex multiplication and complex/real di- 
vision); IH#, OF, UNF set for each result component as appropriate. 





For complex Z, /(Z) is given for the specified functions by the following tables. 


Complex Sine (S1H«Z) 








* 








Argument Z Result 
CHaH “Hah, NaN 3; no exception flags set. 
(t1Inf, *) IVL set; “NaN,NaN> result; message Invalid Arg. 
(*, tInf) RECT Inf, Arg((sin(x), SGH “y2cos(x)))2 »; no exception flags set. 


Ik 





.{F, LHF set for each result component as appropriate. 








Complex Hyperbolic Sine (S 1HH«Z>) 








* 








Argument Z Result 
CHaH “Mah. Hah2; no exception flags set. 
(*, Inf) IVL set; “NaH, NaH? result; message Invalid Ara. 
(+Inf, *) RECT&C Inf, Arg(SGH«x>cosly), sin(y)))® >; no exception flags set. 


IH, OF, LHF set for each result component as appropriate. 














Complex Cosine (i 





Argument Z Result 











CH aH "NaH, Nah; no exception flags set. 
(+Inf, *) IWL set; ¢NaN.HaN> result; message Invalid Arg. 
(*, Inf) RECT &In#. Arg((cos(x), -S “ysin(x)))* >; no exception flags set. 
* IHs, OVF, UHF set for each result component as appropriate. 
] Complex Hyperbolic Cosine (C| €Z>) 
Argument Z Result 





| CHaH “NaN, Mah; no exception flags set. 











(x, Inf) IL set; (NaN, NaH? result; message Invalid Arg. 
(+Inf, *) RECT“ < Inf, Arg((cos(y), 3 «x2 sin(y)))> >; no exception flags set. 
) * IH JF, UNF set for each result component as appropriate. | 





Complex Tangent (TAN <Z)) 





Argument Z Result 
CHaN “Nal. NaN; no exception flags set. 
(EL nt, HINA ca, <y) 3; no exception flags set. 
(tInf, x) IML set; (NaN. NaH > result; message Invalid Ara. 
(*, Inf) ¢ € sin(x)cos(x) 3 #5 , Hy; no exception flags set. 





* O\WF, UNF set for each result component as appropriate. 


Complex Hyperbolic Tangent (TAMH<Z2) 





Argument Z Result 
CH aH “NaN, HaN >; no exception flags set. 
(knit pce Pat) “SGM ¢x> , -; no exception flags set. 
(*, Inf) IML set; “HaH, Nah? result; message Invalid Arg. 
(£0nf, %) : 2; No exception flags set. 





*” Ns, OVF, UNF set for each result component as appropriate. 
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Absolute Value (RES *Z2) 
Argument Z Result 
CHa HH; no exception flags set. 
CiInf Inf; no exception flags set. 
* IH, OMF, UHF set as appropriate. 
Argument (ARG «Z)) 
Argument Z Result 
CHaH Hal; no exception flags set. 
CInf,Inf> 45 degrees or 7/4 radians; IN* set if radian mode. 
C-Int; Inf> 135 degrees or 37/4 radians; IN™* set if radian mode. 
Cinf, ~iné> —45 degrees or —7/4 radians; IH set if radian mode. 
c-Inf.-In#2 | —135 degrees or —3/4 radians; IH® set if radian mode. 


* 


ANGLE*x.y2; INS, UNF set as appropriate. 





Projective Infinity (PFRO1¢Z>) 








* 





Argument Z Result 
CHaH ©MaN, aN; no exception flags set. 
CiInf «Inf, 3; no exception flags set. 


Z; UNF, IH set for any component whose magnitude is between 0 and EF'S if 
TRAPCUHF? # 2. 

















L 











Unit Vector (SGHM¢Z>) 
Argument Z Result 
CHaH "Hah. MsH 2; no exception flags set. 
CZERO Z; no exception flags set. 
(4Inf, Inf) | RECTC“<1, Arg(Z))2; IHX set. 
(+Inf, *) GHex: -y? #69; no exception flags set. 
(x, Inf) -SGH¢x? #8, SGH¢y.> >; no exception flags set. 
* IH, UMF set for each result component as appropriate. 
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Square Root (S@RT«<Z>) 
Argument Z Result 
CH aH “NaH. Ns >; no exception flags set. 
Cinf RECTC¢In#. Arg(Z)“23 3; no exception flags set. 
* IH, UNF set for each result component as appropriate. 








Rectangular to Pola’ 


r Conversion (POLAR ¢Z>) 





Argument Z 


Result 





S 


CABS¢Z),ARG*Z39; INK, OF, UNF set for each result component as 
appropriate. 








Polar to Rectangula’ 


r Conversion (RECT <Z*) 











Argument Z Result 
CHaH “Nah, Nah>; no exception flags set. 
(£in4, & ln) CSGN¢x>+#In¢, 3; no exception flags set. 
(0). SKF) (x, x); no exception flags set. 
(x, Inf) IVL set; “NaH, Nah®> result; message Invalid Ara. 
(+Inf, *) (acos(y),bsin(y)); no exception flags set; 


x if cos(y) + 0 
SGHCx) if cosy) = 0 


=| 


and 
af 


x if sin(y) + 0 
tx. if siny) = 0 


S 








(xcos(y),xsin(y)); IM, LINF set for each result component as appropriate. 











eric Exceptions and the 





Natural Logarithm (LOG«¢Z>) 


T 











Argument Z Result | 
} a 1 
CHat “Hah, NaN; no exception flags set. 
| CZERO et; message LOG (> 
ARGCZ) 
| EAL, ARG A 
CInf (Inf, ARGZ»3; IHX set for the result imaginary part as odienaiie. 


* IH®, UHF set for each result component as appropriate. 


Exponential (Ex 





| Argument Zz | _ Result 














CHall Mah. Nah>; no excesiion ifaw: set. | 
IMaDnt) Bin) &3; no exception flags set. 
(Inf, Inf) | “In#, >; no exception flags set. 
| (x, Inf) IVL set; (NaN.Ns5N% result; message Invalid Arg. 
(-Inf, *) (0 x cos(y),0 x sin(y)); IN set for each result component as appropriate. | 
(Inf, *) RECT£Z2; no exception flags set. 
* | IN#, OWF, LHF set for each vesutt component as appropiate: | 











Relational Operators 


When comparing two walines at least one of which is complex, any numeric comparison operator 
containing « or » without * or # sets I'L and reports message Unordered. If TRAP IWL3 = 2, 
then a result of 0 or 1 will be returned based on the presence of the comparison operator =, that is, Z 
<= W,Z>= W,and Z <>= Ware true if and only if x = uandy = v;Z « W,Z > W,andZ W are 
always false. 





™ (Exponentiation) 


Before W * Z is computed, the following peliminary actions are taken: 


1. If either W or Z is real then, for the purposes of the computation, it becomes complex with 0 
imaginary part. 


nN 


. If either W or Z is a CH ah, then a result of “NaN. HaH> is returned with no exception flags set. 


3. For the purposes of the computation, W and Z are then converted to a canonical form representa- 
tion defined as follows: if one part of a complex number is +Imf while the other part is finite, 
then the canonical form representation replaces the finite part by +0 (that is, preserves its sign); 


otherwise, the complex number is already said to be in canonical form. For example, “4, Inf 3 
and «-Inmf,-—> are the canonical form representations of ©6&,7,I Inf» and ¢-Inf,- 
MAMREAL® respectively. In what follows, W and Z are assumed to be in canonical form. 


For W = CZERO, W ~* Z is given by the table below. 





Exponention (W ~ Z): W = CZERO 


: Argument Z Result 


Cu*x> , 43; no exception flags set. 








x >0 









x <0 set; message B“Heg. 
; E result if TRAF COM 2. 
; C > ZEAL,» result with IH setif TRAP COWS) = 1. 
x = 0 andy = 0 No exception flags set; message 4° reported; default result of «1, 





if TRAPCIVL> #0. | 
O andy # o| IML set; (Ha, NaH? result; message Invalid Ara. | 


x 
I 


For y + 0, W * Z is given by the table below. 


Exponention (W ~ Z): y = 0 


Arguments 
——— Result 
Ww z= 
7 | | = 
j (1, +0) CiInf IL set; “HaH, HaH® result; message i“ Inf. 
* * EXPCZ#LOG¢ (complex x complex multiplication). If Z#LOG°W? 





equals (+0, +Inf), then this quantity is not in the domain of E*P and 
IVL is set, “HaH,HsH> is returned, and message Invalid Ara is 
reported. Otherwise, INS, OWE, and UHF are set for each result compo- 
nent as appropriate. 
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For y = 0 and v # 0, W ~* Z is given by the table below. 


Exponention (W * Z): y = 0 andv+ 0 





Arguments 





IVL set; {HaN,HaN> result; message Invalid Arg. 





CInf CZERO | No exception flags set; message Inf reported; default result of {1,42 
returned if TRAP ¢IWL3 ¥ 0. 
* * EXP¢x#LOG¢W2 > (real + complex multiplication); IN, GVF, UNF set for 


each result component as appropriate. 








For y = 0 and v = 0, W * Z is given by the table below. 


Exponention (W ~ Z): y = 0 andv = 0 





Arguments 





Result 





u=+Inf|x=0 No exception flags set; message Ir“ reported; default result of © 1, 4% 
returned if TRAP <IWL> # 0. 


u= +1 CInf IVL set; (NSaH,MsH> result; message 1“Inf. 
* CInf ©lul*x, 
* * EXP ex#tLOG¢W> ® (real x complex multiplication); I 4% 
each result component as appropriate. 





; No exception flags set. 





IWF, UMF set for 

















Appendix E: Numeric Exceptions and the IEEE Proposal 171 


Array Functions and Operations 


These functions are described in sections 7, 8, and 9 of this manual. Refer to the previous discussion for 
definitions of CZERO, CI nf, complex, etc. 





CHORMCAS, RHORMCA®> 
If Ais M x N (for vectors take N = 1), then 


CHORM 





= . 2 a; RNORMCA) = | MAX = |a;;| 


patAX =j<=M 


If any element of A is a signaling Hal (inchiding either part of complex AEA elements), then each 
function sets I'L and reports message Signaled Op. If TRAPS ITVL} = 2, the result is a quiet 
‘ah with no other elements processed. 


If any element of A is a quiet 14N (including either part of complex array elements), then each func- 
tion sets IL and reports message Unior dered; a Hat result is returned. Otherwise, IM, O\F, and 
LHF are set for the result as appropriate. 


FHORM* A? 


If Ais M x N (for vectors take N = 1), then 


FHORMCA) = 





If any element of A is a signaling !/34 (including either part of complex array elements), then J. is 
set and message Signaled Op is reported. If TRAP“ I\'L% = 2, the result is a quiet 4a with no 
other elements processed. 


Quiet als pass through with no exception flags set. Otherwise, 1M, O\/F, and LIF are set for the 
result as appropriate. 





OOTCA, Ba 


If A and B are N-element vectors, then 

N -_. 

DOT<A.B) = & ap; 

i=1 
(If either A or B is complex, refer to the definitions of complex addition and multiplication given 
previously). If any element of A or B is a signaling 43H (including either part of complex array ele- 
ments), then I'L is set with message Signaled Op. If, in any term in the above expression, +0 or 
CZERO is multiplied by +Inf or CIn#, then I\/L is set with message Inf +4. If, in the above 
expression, the summation executes an addition equivalent to Inf — Inf, then I\'L is set with 
message Inf-Inf. 


If only one I\/L exception occurs, that message is reported. If more than one I\'L exception occurs, 
the particular message(s) reported depends upon the order and type of exception that occurs. If 
TRAP CIWL3 = 2, the result is either a real Hat or a complex value with one or two Hah compo- 
nents. Quiet ahs pass through with no exception flags set. Otherwise, IN, OF, and UMF are set for 
the result, or each result component, as appropriate. 


MAT C=AB 
If Ais M x Nand Bis N x P (for vectors take P=1), then 
N 
cy = = indy; 
k=1 


(If either A or B is complex, refer to the definitions of complex addition and multiplication given 
previously). Since each result element is derived from an inner product, exception handling is the same 
as that for 00 T< A.B», applied to each result element separately. 


MAT C=TEH CA #B 


If Ais M x N and Bis M x P (for vectors take P=1), then 


M 


ci = 2 pid; 


(If either A or B is complex, refer to the definitions of complex addition and multiplication given 
previously). 


Since each result element is derived from an inner product, exception handling is the same as that for 
O0T*¢A,B2, applied to each result element separately. 
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MAT C=A+B 


All elements of C are computed separately as 


(If either A or B is complex, refer to the definitions of complex addition and subtraction given pre- 
viously). 

If any element of A or B is a signaling HH (including either part of complex array elements), then 
IVL is set and message Signaled Op is reported. If TRAF «IL» = 2, the corresponding result 
element or component becomes a quiet 13M and the operation continues. Quiet als pass through 
with no exception flags set. 

IL is set and message Im#~Inrf is reported if any addition or subtraction (or componentwise addi- 


tion or subtraction) is equivalent to Inf — Inf; a NaN is returned for the corresponding result 
element or component. Otherwise, IN, O\/F, and UNF are set for each result element or component 


as appropriate. 


NAT B=tsttA 
All elements of B are computed separately as 
by = say 


(If either s or A is complex, refer to the definition of complex multiplication given previously). If s (or 
either part of s, if s is complex) is a signaling NaN, then I\/L is set and message Signaled Op is 
eported; if any element of A is a signaling HN (including either part of complex array elements), then 
IL is set and message Signaled (ip is reported. In either event, if TRAP ¢IVL» = 2, these 
Hats become quiet and the operation continues. Quiet als pass through with no exception flags set. 


IL is set and message If #4 is reported if, during the computation of any result element, +0 or 
CZERO is multiplied by + Inf or Clnf. If TRAP’ IL = 2, the corresponding result element is 
either a real HS or a complex value with one or two Hs components. Otherwise, IMs, OF, and 
UNF are set for each result element or component as appropriate. 


DETCA?, MAT C=IHVCA>, HAT C=5 [A.B 
Due to the intricate algorithmic basis of these three operations, exception handling is complex; only a 
ummary is provided here. 








any element of A or B is a signaling =H (including either part of complex array elements), then 
VL is set and message Signaled Op is reported. If TRAP“ I'L 3 = 2, the corresponding element 
component becomes a quiet Ha and the operation continues. Quiet Has pass through with no 


kception flags set. 
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OWF, UNF, and IN® are set for each result element as appropriate and may also be set at intermediate 
stages of the computation (especially 0'/F when A is (machine) singular). I‘! may also be set with 
any of the following messages reported: Inf #9, Imf—In#, and/or Inf«Inf. These messages are 
only possible due to a + Inf in A or B or an intermediate overflow becoming + If; in the latter case 
they may be suppressed by setting TRAP «<G'F > = 1 before the computation. 


Other Math Pac Functions 
PROOT 


Special cases for the FF OCT function are handled first. These are }}3's, Infs, or leading and trailing 
zeros in the coefficient array. 


Nats are handled first. If any coefficient is a Hah, then every element of the result array becomes 
“Mah. Hah with no exception flags set and the function is complete. (Signaling aN coefficients do 
not set T/L). 


Ins are dealt with next. If any coefficient is + If, then every finite coefficient will become zero and 
the computation falls through to handle leading and trailing zeros. 


Leading zeros are handled next. Every leading zero coefficient will produce a root at “Inf, Inf 3 
with no exception flags set. The next coefficient then becomes the leading coefficient and the process 
loops. Every such root stored decrements the degree of the polynomial; the function is complete if the 
degree becomes zero. 


Trailing zeros are handled next. Every trailing zero coefficient will produce a root at ‘4, &% with no 
exception flags set. The second to the last coefficient then becomes the trailing coefficient and the 
process loops. Every such root stored decrements the degree of the polynomial and the function is 
complete if the degree becomes zero. 


At this point, the degree of the polynomial is positive and either all (remaining) coefficients are finite, 
in which case the roots of the (reduced) polynimial will be found, or the leading and trailing coefficients 
are both +Ir#. In the latter case, at least two of the original coefficients were +11f and factoriza- 
tion does not make sense; if the (new) degree of the polynomial is D, then D roots at {HaH,.MaH? are 
stored into the result array and the function is complete; every such root stored sets I'L and reports 
message Invalid Ara. 


Except for the above special cases, 'F and LHF are set for every result array component as appro- 
priate with IH always set. 
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FOUR 


As with the FET function, special cases for the FOF function are handled first. These are Hal 
and Inf components in the data array. 


‘ahs are handled first. If any component of any data array element is a Mati, then every element of 
the result array becomes “Ha, HaH> with no exception flags set and the function is complete. 
(Signaling 44H components do not set I'L). 


Imfs are dealt with next. If any component of any data array element is +Imf, then every result 
element becomes «In? . Inf? with no exception flags set and the function is complete. 


Except for the above special cases, 0'!F and LIHF are set for every result array component as appro- 
priate with I! always set unless the data array was identically zero. 


FHREOOT and THTEGRAL 


If a aN (signaling or quiet) results during the evaluation of any of the arguments of FMR OOT or 
IHTEGRAL, then error Invalid Arg is reported; no exception flags are set and this error halts 
the computation. 


In general, any value of +I resulting from the evaluation of any of the arguments of FHROOT or 
IHTEGEAL becomes +MASREAL for the purposes of the computation. IH, O\F, and LMF are set 
for the result as appropriate. 


Remember that FMPOCT looks at the value of TRAP «UNF» to decide whether or not to search the 
range of denormalized numbers for a root. This region is searched only if TRAF ¢\JMF » = 2 when the 
FHROOT function is started. 





Keyword Index 


= oe = —— = 


| Page | Description 


41‘ Absolute value of a complex number. 

28 Inverse hyperbolic cosine. 

41 Argument of a complex number. 

28 Inverse hyperbolic sine. 

28 Inverse hyperbolic tangent. 

16 Decimal to binary/octal/hexadecimal conversion. 
15 _—_ Binary/octal/hexadecimal to decimal conversion. 
es Srey 22 Complex IMAGE field. 

CHORM 70 One-norm (column norm) of an array. 
COMPLEX 19 Complex variable creation. 

SHORT 19 Complex short variable creation. 

21 Conversion, real to complex. 

42 Complex conjugate. 

38 Complex cosine. 

27 Hyperbolic cosine. 

39 Complex hyperbolic cosine. 

69 Determinant of a matrix. 





xn 
imo 

x 
we 











CET (no operand) 69 Determinant of last real matrix used as operand of IN|! or first 
operand of S15. 
DETL 69 Same as LIET (no operand). 


ooT 71 Dot (inner) product. 
: 37 Complex exponential (e7) 





} 90  Second-best guess to value returned by last FMR OMT. 
FHORM 70 ~~‘ Frobenius norm. 





FHROOT 89  Rootfinding for functions. | 

FMALUE 90 Functional value of last FHROOT. | 
| FUAR 90  ~=Variable to solve for in FHROOT. 
| GAMMA 28 Gamma function. 

IBOUNO 103. -Uncertainty of last INTEGRAL. 

IMFT 21 Imaginary part of complex number. 

INTEGRAL 101 Integration of functions. 

TROUHO 30 ~~ Integer round. 
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Keyword 
IVALUE 
| IVAR 

LENO 
| LBQUHO 
| Loc 
Loge 

MAT OISP 
MAT O1 
| MAT INPLIT 
MAT «1. COHN 
| MAT... TDN 
) MAT... ZER 
MAT,..2ERO 
| MAT. ..PRINT 
| (MAT 
| MAT = 
| MAT =- 
MAT a9 
MAT =...- 
MAT =. # 
MAT =¢3 
| MAT =co# 
HAT = FOUR 
MAT = IMM 
MAT = PROOT 
MAT 
| MAT = 
| MAT = 
| HAW 
HETGHBOR 












USTHG 








TRH 
TERA a F 














PRINT USIHG 


135 
77 
120 
73. 
LT 
66 
30 
30 
40 
42 
40 
21 
70 
29 
41 
38 


Description 


Current approximation to an IHTE 
Variable of integration in INTEGRAL. 
Array subscript lower bound. 

Same as LENO. 

Complex natural logarithm. 

Log base 2. 

Array display (unformatted). 

Array display (formatted). 

Interactive array input. 

Constant array with redimensioning. 
Identity matrix with redimensioning. 
Zero array with redimensioning. 

Same as MAT... ZER. 

Array printing (unformatted). 

Array printing (formatted). 

Array copying (simple assignment). 
Array negation. 

Array addition. 

Array subtraction. 

Array multiplication. 





Scalar to array assignment (numeric expression assignment). 


Scalar multiplication. 

Finite Fourier Transform. 

Matrix inversion. 

Polynomial rootfinding. 

System solution. 

Transpose or conjugate transpose. 
Transpose or conjugate transpose multiply. 
14H diagnostic function. 
Successor/predecessor function. 
Rectangular to polar conversion. 
Conversion of complex infinities to projective infinities. 
Polar to rectangular conversion. 

Real part of complex number. 

Infinity (row) norm of an array. 

Exponent scaling function. 

Complex unit vector. 

Complex sine. 
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Keyword | Page Description 
27 Hyperbolic sine. 
39 Complex hyperbolic sine. 
40 Complex square root. 
40 Same as SOR. 
38 Complex tangent. 
TAHH 27 Hyperbolic tangent. 
TAHH 39 Complex hyperbolic tangent. 
TYPE 31 Data type function. 
WEHO 71. Array subscript upper bound. 
UBOUND 71 Same as WEHE. 
+ 35 Complex addition. 
= 35 Complex unary minus. 
= 36 Complex subtraction. 
4 36 Complex multiplication. 
36 Complex division. 
36 Complex exponentiation (Z”) 
y 43 Complex relational operators. 
# 
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